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Abstract 

Recent cosmological observations suggest the existence of a positive cosmological 
constant A with the mag nitude K{Gh/c') ^ IQ-^^a^ r^^A^ 

review discusses several 

aspects of the cosmological constant both from the cosmological (sections 1-6) and 
field theoretical (sections 7-11) perspectives. After a brief introduction to the key 
issues related to cosmological constant and a historical overview, a summary of 
the kinematics and dynamics of the standard Friedmann model of the universe is 
provided. The observational evidence for cosmological constant, especially from the 
supernova results, and the constraints from the age of the universe, structure for- 
mation, Cosmic Microwave Background Radiation (CMBR) anisotropics and a few 
others are described in detail, followed by a discussion of the theoretical models 
(quintessence, tachyonic scalar field, ...) from different perspectives. The latter part 
of the review (sections 7-11) concentrates on more conceptual and fundamental 
aspects of the cosmological constant like some alternative interpretations of the 
cosmological constant, relaxation mechanisms to reduce the cosmological constant 
to the currently observed value, the geometrical structure of the de Sitter space- 
time, thermodynamics of the de Sitter universe and the role of string theory in the 
cosmological constant problem. 
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1 Introduction 

This review discusses several aspects of the cosmological constant both from 
the cosmological and field theoretical perspectives with the emphasis on con- 
ceptual and fundamental issues rather than on observational details. The plan 
of the review is as follows: This section introduces the key issues related to cos- 
mological constant and provides a brief historical overview. (For previous re- 
views of this subject, from cosmological point of view, see [1,2,3,139].) Section 

2 summarizes the kinematics and dynamics of the standard Friedmann model 
of the universe paying special attention to features involving the cosmological 
constant. Section 3 reviews the observational evidence for cosmological con- 
stant, especially the supernova results, constraints from the age of the universe 
and a few others. We next study models with evolving cosmological 'constant' 
from different perspectives. (In this review, we shall use the term cosmological 
constant in a generalized sense including the scenarios in which cosmological 
"constant" is actually varying in time.) A phenomenological parameterization 
is introduced in section 4.1 to compare theory with observation and is followed 
up with explicit models involving scalar fields in section 4.2. The emphasis is 
on quintessence and tachyonic scalar field models and the cosmic degeneracies 
introduced by them. Section 5 discusses cosmological constant and dark en- 
ergy in the context of models for structure formation and section 6 describes 
the constraints arising from CMBR anisotropies. 

The latter part of the review concentrates on more conceptual and fundamen- 
tal aspects of the cosmological constant. ( For previous reviews of this subject, 
from a theoretical physics perspective, see [4,5,6].) Section 7 provides some 
alternative interpretations of the cosmological constant which might have a 
bearing on the possible solution to the problem. Several relaxation mechanisms 
have been suggested in the literature to reduce the cosmological constant to 
the currently observed value and some of these attempts are described in sec- 
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tion 8. Section 9 gives a brief description of the geometrical structure of the 
de Sitter spacetime and the thermodynamics of the de Sitter universe is taken 
up in section 10. The relation between horizons, temperature and entropy are 
presented at one go in this section and the last section deals with the role of 
string theory in the cosmological constant problem. 



1.1 The many faces of the cosmological constant 



Einstein's equations, which determine the dynamics of the spacetime, can be 
derived from the action (see, eg. [7]): 

A = I RV^d^'x + I L^atter(0, d(t>)yf^d''x (l) 



where Lmatter is the Lagrangian for matter depending on some dynamical vari- 
ables generically denoted as 0. (We are using units with c—1.) The variation 
of this action with respect to will lead to the equation of motion for matter 
(5Linatter/^0) = 0, in a given background geometry, while the variation of the 
action with respect to the metric tensor leads to the Einstein's equation 

Rik - \9^kR = levrG^^^ = STrGTife (2) 



where the last equation defines the energy momentum tensor of matter to be 

T,fc=2(5L^,tter/5<7*')- 

Let us now consider a new matter action -^^^attcr ~ -^matter ~ (A/SttG) where A 
is a real constant. Equation of motion for the matter (^Lmattcr/^^) = 0, does 
not change under this transformation since A is a constant; but the action 
now picks up an extra term proportional to A 



^ = 16^ / + / - 8^) ^ 



(3) 



and equation (2) gets modified. This innocuous looking addition of a constant 
to the matter Lagrangian leads to one of the most fundamental and fascinating 
problems of theoretical physics. The nature of this problem and its theoretical 
backdrop acquires different shades of meaning depending which of the two 
forms of equations in (3) is used. 

The first interpretation, based on the first line of equation (3), treats A as the 
shift in the matter Lagrangian which, in turn, will lead to a shift in the matter 



4 



Hamiltonian. This could be thoTight of as a shift in the zero point energy of 
the matter system. Such a constant shift in the energy does not affect the 
dynamics of matter while gravity — which couples to the total energy of the 
system — picks up an extra contribution in the form of a new term Qjfc in the 
energy-momentum tensor, leading to: 

Rl - l^iR = 87rG(T^ + Ql); Q\ ^ ^ (4) 

The second line in equation (3) can be interpreted as gravitational field, de- 
scribed by the Lagrangian of the form Lgrav oc — 2A), interacting with 
matter described by the Lagrangian ^matter- In this interpretation, gravity is 
described by two constants, the Newton's constant G and the cosmological 
constant A. It is then natural to modify the left hand side of Einstein's equa- 
tions and write (4) as: 

Rl - I^IR - SIA = MGTl (5) 

In this interpretation, the spacetime is treated as curved even in the absence of 
matter {Tik = 0) since the equation Rik — {l/2)gikR — Agik = does not admit 
flat spacetime as a solution. (This situation is rather unusual and is related 
to the fact that symmetries of the theory with and without a cosmological 
constant are drastically different; the original symmetry of general covariance 
cannot be naturally broken in such a way as to preserve the sub group of 
spacetime translations.) 

In fact, it is possible to consider a situation in which both effects can occur. 
If the gravitational interaction is actually described by the Lagrangian of the 
form (i?— 2A), then there is an intrinsic cosmological constant in nature just as 
there is a Newtonian gravitational constant in nature. If the matter Lagrangian 
contains energy densities which change due to dynamics, then Lmatter can pick 
up constant shifts during dynamical evolution. For example, consider a scalar 
field with the Lagrangian -Lmatter = (l/2)9i09V — which has the energy 
momentum tensor 

T« = - 6t - (6) 

For field configurations which are constant [occurring, for example, at the 
minima of the potential 1^(0)], this contributes an energy momentum tensor 
'^b ~ ^6^(0min) which has exactly the same form as a cosmological constant. 
As far as gravity is concerned, it is the combination of these two effects — 
of very different nature — which is relevant and the source will be T^f = 
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[l^(0min) + (A/87rG)]g'a&, corresponding to an effective gravitational constant 



If 0min and lience V(0min) clianges during dynamical evolution, the value of 
Aeff can also change in course of time. More generally, any field configuration 
which is varying slowly in time will lead to a slowly varying Aeg. 

The extra term Qik in Einstein's equation behaves in a manner which is very 
peculiar compared to the energy momentum tensor of normal matter. The 
term Ql. = paS^ is in the form of the energy momentum tensor of an ideal 
fluid with energy density pA and pressure Pa — —pA] obviously, either the 
pressure or the energy density of this "fluid" must be negative, which is unlike 
conventional laboratory systems. (See, however, reference [8].) 

Such an equation of state, p = —P also has another important implication 
in general relativity. The spatial part g of the geodesic acceleration (which 
measures the relative acceleration of two geodesies in the spacctimc) satisfies 
the following exact equation in general relativity (see e.g., page 332 of [9]): 



showing that the source of geodesic acceleration is {p + 3P) and not p. As 
long as (p + 3P) > 0, gravity remains attractive while {p + 3P) < can 
lead to repulsive gravitational effects. Since the cosmological constant has 
(Pa + 3Pa) = —2pA, a positive cosmological constant (with A > 0) can lead to 
repulsive gravity. For example, if the energy density of normal, non-relativistic 
matter with zero pressure is pnr, then equation (8) shows that the geodesies 
will accelerate away from each other due to the repulsion of cosmological 
constant when pnr < 2pA- A related feature, which makes the above con- 
clusion practically relevant is the fact that, in an expanding universe, pA 
remains constant while pnr decreases. (More formally, the equation of mo- 
tion, d{pAV) = —PAdV for the cosmological constant, treated as an ideal 
fluid, is identically satisfied with constant Pa,Pa-) Therefore, pa will eventu- 
ally dominate over pnr if the universe expands sufficiently. Since |Ap/^ has 
the dimensions of inverse length, it will set the scale for the universe when 
cosmological constant dominates. 

It follows that the most stringent bounds on A will arise from cosmology when 
the expansion of the universe has diluted the matter energy density sufficiently. 

The rate of expansion of the universe today is usually expressed in terms of 
the Hubble constant: Hq = lOOh km s^^ Mpc^^ where 1 Mpc ^ 3 x 10^^ cm 
and /i is a dimensionless parameter in the range 0.62 < /i < 0.82 (see section 
3.2). Prom Hq we can form the time scale t^niv = Hq^ 10^°/i~^ yr and 



Aeff = A + 87rGV(0^in) 



(7) 



V • g = -inGip + 3P) 



(8) 
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the length scale cHq^ ^ 3000/i~^ Mpc; tuniv characterizes the evolutionary 
time scale of the universe and H^^ is of the order of the largest length scales 
currently accessible in cosmological observations. From the observation that 
the universe is at least of the size Hq^, we can set a bound on A to be 
|A| < 10"^^ cm~^. This stringent bound leads to several issues which have 
been debated for decades without satisfactory solutions. 

• In classical general relativity, based on the constants G, c and A, it is not 
possible to construct any dimensionless combination from these constants. 
Nevertheless, it is clear that A is extremely tiny compared to any other phys- 
ical scale in the universe, suggesting that A is probably zero. We, however, 
do not know of any symmetry mechanism or invariance principle which re- 
quires A to vanish. Supersymmetry does require the vanishing of the ground 
state energy; however, supersymmetry is so badly broken in nature that this 
is not of any practical use [10,11]. 

• We mentioned above that observations actually constrain Acr in equation 
(7), rather than A. This requires A and V(0mm) to be fine tuned to an 
enormous accuracy for the bound, |Aeff| < 10~®^ cm~^, to be satisfied. This 
becomes more mysterious when we realize that V(0min) itself could change 
by several orders of magnitude during the evolution of the universe. 

• When quantum fields in a given curved spacetime are considered (even with- 
out invoking any quantum gravitational effects) one introduces the Planck 
constant, h, in the description of the physical system. It is then possible to 
form the dimensionless combination KiGh/c") = ALp. (This equation also 
defines the quantity Lp; throughout the review we use the symbol '=' to de- 
fine variables.) The bound on A translates into the condition ALp < lO""*^^^. 
As has been mentioned several times in literature, this will require enormous 
fine tuning. 

• All the above questions could have been satisfactorily answered if we take 
Aeff to be zero and assume that the correct theory of quantum gravity will 
provide an explanation for the vanishing of cosmological constant. Such a 
view was held by several people (including the author) until very recently. 
Current cosmological observations however suggests that Acfr is actually 
non zero and AggLp is indeed of order O{10^^^^). In some sense, this is the 
cosmologist's worst nightmare come true. If the observations are correct, 
then Aeff is non zero, very tiny and its value is extremely fine tuned for no 
good reason. This is a concrete statement of the first of the two 'cosmological 
constant problems'. 

• The bound on ALp arises from the expansion rate of the universe or — 
equivalently — from the energy density which is present in the universe to- 
day. The observations require the energy density of normal, non relativistic 
matter to be of the same order of magnitude as the energy density con- 
tributed by the cosmological constant. But in the past, when the universe 
was smaller, the energy density of normal matter would have been higher 
while the energy density of cosmological constant does not change. Hence 



7 



we need to adjust the energy densities of normal matter and cosmological 
constant in the early epoch very carefully so that pA ^ Pnr around the 
current epoch. If this had happened very early in the evolution of the uni- 
verse, then the repulsive nature of a positive cosmological constant would 
have initiated a rapid expansion of the universe, preventing the formation of 
galaxies, stars etc. If the epoch of p\ ~ pnr occurs much later in the future, 
then the current observations would not have revealed the presence of non 
zero cosmological constant. This raises the second of the two cosmological 
constant problems: Why is it that (pa/pnr) = at the curreni phase of 
the universe ? 

• The sign of A determines the nature of solutions to Einstein's equations as 
well as the sign of {p\ + 3Pa). Hence the spacetime geometry with ALp = 
"1^0-123 jg YQTj-y different from the one with ALp = — 10~^^^. Any theoretical 
principle which explains the near zero value of ALp must also explain why 
the observed value of A is positive. 

At present we have no clue as to what the above questions mean and how 
they need to be addressed. This review summarizes different attempts to un- 
derstand the above questions from various perspectives. 



1.2 A brief history of cosmological constant 

Originally, Einstein introduced the cosmological constant A in the field equa- 
tion for gravity (as in equation (5)) with the motivation that it allows for a 
finite, closed, static universe in which the energy density of matter determines 
the geometry. The spatial sections of such a universe are closed 3-spheres with 
radius / = (SttG^nr)^^^^ = A~^/^ where pnr is the energy density of pres- 
sureless matter (see section 2.4) Einstein had hoped that normal matter is 
needed to curve the geometry; a demand, which — to him — was closely re- 
lated to the Mach's principle. This hope, however, was soon shattered when de 
Sitter produced a solution to Einstein's equations with cosmological constant 
containing no matter [12]. However, in spite of two fundamental papers by 
Priedmann and one by Lemaitre [13,14], most workers did not catch on with 
the idea of an expanding universe. In fact, Einstein originally thought Pried- 
mann's work was in error but later published a retraction of his comment; 
similarly, in the Solvay meeting in 1927, Einstein was arguing against the so- 
lutions describing expanding universe. Nevertheless, the Einstein archives do 
contain a postcard from Einstein to Weyl in 1923 in which he says: "If there is 
no quasi-static world, then away with the cosmological term" . The early his- 
tory following de Sitter's discovery (see, for example, [15]) is clearly somewhat 
confused, to say the least. 

It appears that the community accepted the concept of an expanding universe 
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largely due to the work of Lemaitre. By 1931, Einstein himself had rejected the 
cosmological term as superflous and imjustificd (see reference [16], which is a 
single authored paper; this paper has been mis-cited in literature often, even- 
tually converting part of the journal name "preuss" to a co-author "Preuss, S. 
B"!; see [17]). There is no direct record that Einstein ever called cosmological 
constant his biggest blunder. It is possible that this often repeated "quote" 
arises from Gamow's recollection [18]: "When I was discussing cosmological 
problems with Einstein, he remarked that the introduction of the cosmological 
term was the biggest blunder he ever made in his life." By 1950's the view 
was decidedly against A and the authors of several classic texts ( like Landau 
and Liftshitz[7], Pauh [19] and Einstein [20]) argued against the cosmological 
constant. 

In later years, cosmological constant had a chequered history and was often 
accepted or rejected for wrong or insufficient reasons. For example, the original 
value of the Hubble constant was nearly an order of magnitude higher [21] than 
the currently accepted value thereby reducing the age of the universe by a 
similar factor. At this stage, as well as on several later occasions (eg., [22,23]), 
cosmologists have invoked cosmological constant to reconcile the age of the 
universe with observations (see section 3.2). Similar attempts have been made 
in the past when it was felt that counts of quasars peak at a given phase in 
the expansion of the universe [24,25,26]. These reasons, for the introduction 
of something as fundamental as cosmological constant, seem inadequate at 
present. 

However, these attempts clearly showed that sensible cosmology can only be 
obtained if the energy density contributed by cosmological constant is compa- 
rable to the energy density of matter at the present epoch. This remarkable 
property was probably noticed first by Bondi [27] and has been discussed by 
McCrea [28]. It is mentioned in [1] that such coincidences were discussed in 
Dicke's gravity research group in the sixties; it is almost certain that this must 
have been noticed by several other workers in the subject. 

The first cosmological model to make central use of the cosmological con- 
stantwas the steady state model [29,30,31]. It made use of the fact that a 
universe with a cosmological constant has a time translational invariance in a 
particular coordinate system. The model also used a scalar field with negative 
energy field to continuously create matter while maintaining energy conser- 
vation. While modern approaches to cosmology invokes negative energies or 
pressure without hesitation, steady state cosmology was discarded by most 
workers after the discovery of CMBR. 

The discussion so far has been purely classical. The introduction of quantum 
theory adds a new dimension to this problem. Much of the early work [32,33] 
as well as the definitive work by Pauli [34,35] involved evaluating the sum 
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of the zero point energies of a quantum field (with some cut-off) in order 
to estimate the vacuum contribution to the cosmological constant. Such an 
argument, however, is hopelessly naive (inspite of the fact that it is often 
repeated even today). In fact, Pauli himself was aware of the fact that one must 
exclude the zero point contribution from such a calculation. The first paper to 
stress this clearly and carry out a second order calculation was probably the 
one by Zeldovich [36] though the connection between vacuum energy density 
and cosmological constant had been noted earlier by Gliner [37] and even 
by Lemaitre [38]. Zeldovich assumed that the lowest order zero point energy 
should be subtracted out in quantum field theory and went on to compute 
the gravitational force between particles in the vacuum fluctuations. If E is 
an energy scale of a virtual process corresponding to a length scale / = he/ 
then — [E/hcY particles per unit volume of energy E will lead to the 
gravitational self energy density of the order of 

G{E/c^Y^ 3 GE^ 

This wiU correspond to AL|, ^ {E/Epf where Ep = {hc^/Gf^ ^ lO^^GeV 
is the Planck energy. Zeldovich took E Ri 1 GeV (without any clear reason) 
and obtained a pa which contradicted the observational bound "only" by nine 
orders of magnitude. 

The first serious symmetry principle which had implications for cosmologi- 
cal constant was supersymmetry and it was realized early on [10,11] that the 
contributions to vacuum energy from fermions and bosons will cancel in a su- 
persymmetric theory. This, however, is not of much help since supersymmetry 
is badly broken in nature at sufficiently high energies (at Eqs > 10^ Gev). 
In general, one would expect the vacuum energy density to be comparable to 
the that corresponding to the supersymmetry braking scale. Ess- This will, 
again, lead to an unacceptably large value for p\. In fact the situation is more 
complex and one has to take into account the coupling of matter sector and 
gravitation — which invariably leads to a supergravity theory. The description 
of cosmological constant in such models is more complex, though none of the 
attempts have provided a clear direction of attack (see e.g, [4] for a review of 
early attempts). 

The situation becomes more complicated when the quantum field theory ad- 
mits more than one ground state or even more than one local minima for the 
potentials. For example, the spontaneous symmetry breaking in the electro- 
weak theory arises from a potential of the form 

y = yo-/xV + #' (/^',y>o) (10) 

At the minimum, this leads to an energy density V^m = Vq— (/i^/45f). If we take 
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Vo = then {Vrmn/ g) ~ —(300 GeV)^. For an estimate, we will assume that 
the gauge coupling constant g is comparable to the electromagnetic coupling 
constant: g = 0{a'^), where a = (e^/hc) is the fine structure constant. Then, 
we get \Vrnin\ ~ 10^ GeV^ which misses the bound on A by a factor of 10^^. It is 
really of no help to set Vmin = by hand. At early epochs of the universe, the 
temperature dependent effective potential [39,40] will change minimum to (f) — 
with V{(j)) = Vq. In other words, the ground state energy changes by several 
orders of magnitude during the electro-weak and other phase transitions. 

Another facet is added to the discussion by the currently popular models 
of quantum gravity based on string theory [41,42]. The currently accepted 
paradigm of string theory encompasses several ground states of the same un- 
derlying theory (in a manner which is as yet unknown). This could lead to 
the possibility that the final theory of quantum gravity might allow different 
ground states for nature and we may need an extra prescription to choose the 
actual state in which we live in. The different ground states can also have differ- 
ent values for cosmological constant and we need to invoke a separate (again, 
as yet unknown) principle to choose the ground state in which ALp 10~^^^ 
(see section 11). 



2 Framework of standard cosmology 



All the well developed models of standard cosmology start with two basic 
assumptions: (i) The distribution of matter in the universe is homogeneous 
and isotropic at sufficiently large scales, (ii) The large scale structure of the 
universe is essentially determined by gravitational interactions and hence can 
be described by Einstein's theory of gravity. The geometry of the universe 
can then be determined via Einstein's equations with the stress tensor of 
matter T^(t,x) acting as the source. (For a review of cosmology, see e.g. 
[43,44,45,46,47]). The first assumption determines the kinematics of the uni- 
verse while the second one determines the dynamics. We shall discuss the 
consequences of these two assumptions in the next two subsections. 



2. 1 Kinematics of the Friedmann model 



The assumption of isotropy and homogeneity implies that the large scale ge- 
ometry can be described by a metric of the form 



dt^ - a\t)dyi'' ^ dt^ - a\t) 



dr^ 
1 — kr"^ 
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in a suitable set of coordinates called comoving coordinates. Here a{t) is an 
arbitrary function of time (called expansion factor) and k = 0, ±1. Defining 
a new coordinate x through x = (r, sin~^ r, sinh~^ r) for k = (0, +1, —1) this 
line element becomes 



els'" = df - a'^dx' = df - a'(i) [dx^ + SHx) (dO^ + sin^ ■ 



ed(j)'^)\ (12) 



where Sk{x) — (x, sinx, sinhx) for k — (0, +1,— 1). In any range of time 
during which a{t) is a monotonic function of t, one can use a itself as a time 
coordinate. It is also convenient to define a quantity called the redshift, 
through the relation a{t) = ao[l + z{t)]~^ where oq is the current value of the 
expansion factor. The line element in terms of [a, x, 0,4'] or %, ^, 0] is: 



ds^ = H-\a) 



H-\z)dz' - dx' 



(13) 



where H{a) = (a/a), called the Hubble parameter, measures the rate of ex- 
pansion of the universe. 

This equation allows us to draw an important conclusion: The only non trivial 
metric function in a Friedmann universe is the function H{a) (and the numer- 
ical value of k which is encoded in the spatial part of the line element.) Hence, 
any kind of observation based on geometry of the spacetime, however com- 
plex it may be, will not allow us to determine anything other than this single 
function H{a). As we shall see, this alone is inadequate to describe the mate- 
rial content of the universe and any attempt to do so will require additional 
inputs. 

Since the geometrical observations often rely on photons received from distant 
sources, let us consider a photon traveling a distance rem (-2) from the time 
of emission (corresponding to the redshift z) till today. Using the fact that 
ds = for a light ray and the second equality in equation (13) we find that 
the distance traveled by light rays is related to the redshift by dx = H~^{z)dz. 
Integrating this relation, we get 

1 r 

rem(^) = Skia); a=— H-\z)dz (14) 

ao^ 



All other geometrical distances can be expressed in terms of r^^^z) (see eg., 
[44]). For example, the flux of radiation F received from a source of luminosity 
L can be expressed in the form F — L/ [A.T:d\) where 

dhiz) = aorem(^)(l + z) = ao(l + z)Sk{oL) (15) 
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is called the luminosity distance. Similarly, if D is the physical size of an object 
which subtends an angle 6 to the observer, then — for smah 5 — we can define 
an angular diameter distance dA through the relation 5 = D/dA- The angular 
diameter distance is given by 

dA{z) = aorem(^)(l + z)~^ (16) 



with dL — (l + zYdA- 

If we can identify some objects (or physical phenomena) at a redshift of z 
having a characteristic transverse size £), then measuring the angle 5 sub- 
tended by this object we can determine dA{z). Similarly, if wc have a series 
of luminous sources at different redshifts having known luminosity L, then 
by observing the flux from these sources L, one can determine the luminosity 
distance dL^z). Determining any of these functions will allow us to use the 
relations (15) [or (16)] and (14) to obtain H~^{z). For example, H~^{z) is 
related to diiz) through 



H-\z) 



1 - 



kdl{z) 


d 


diiz) 


d 


diiz) 


al{l + zf\ 


dz 


.1 + -^. 


dz 


.1 + -^. 



(17) 



where second equality holds if the spatial sections of the universe are flat, 
corresponding to A; = 0; then dL{z), dA{z),rem{z) and H~^{z) all contain the 
(same) maximal amount of information about the geometry. 



The function Vcmiz) also determines the proper volume of the universe between 
the redshifts z and z+dz subtending a solid angle dfl in the sky. The comoving 
volume element can be expressed in the form 



dV 

didn ^ ""^^"^ 



dr 
dz 



dl 



oc 



[i + zY 



[I + z)d'L 

dj. 



- 1 



(18) 



where the prime denotes derivative with respect to z. Based on this, there 
has been a suggestion [48] that future observations of the number of dark 
matter halos as a function of redshift and circular velocities can be used to 
determine the comoving volume element to within a few percent accuracy. If 
this becomes possible, then it will provide an additional handle on constraining 
the cosmological parameters. 

The above discussion illustrates how cosmological observations can be used to 
determine the metric of the spacetime, encoded by the single function H~^{z). 
This issue is trivial in principle, though enormously complicated in practice 
because of observational uncertainties in the determination of (ix,(2;), dA^z) etc. 
We shall occasion to discuss these features in detail later on. 
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2.2 Dynamics of the Friedmann model 

Let us now turn to the second assumption which determines the dynamics of 
the universe. When several non interacting sources are present in the universe, 
the total energy momentum tensor which appear on the right hand side of 
the Einstein's equation will be the sum of the energy momentum tensor for 
each of these sources. Spatial homogeneity and isotropy imply that each 
is diagonal and has the form — dia [pi(t), — Pj(t), — Pj(t), — Pj(t)] where 
the index i — 1,2, ...,N denotes different kinds of sources (like radiation, 
matter, cosmological constant etc.). Since the sources do not interact with 
each other, each energy momentum tensor must satisfy the relation T^"^ = 
which translates to the condition d^piO?) = —Pida^. It follows that the 
evolution of the energy densities of each component is essentially dependent 
on the parameter Wi = {Pi/ pi) which, in general, could be a function of time. 
Integrating d{pia^) — —Wipida^, we get 




which determines the evolution of the energy density of each of the species in 
terms of the functions Wi{a). 

This description determines p{a) for different sources but not a{t). To deter- 
mine the latter we can use one of the Einstein's equations: 

i 

This equation shows that, once the evolution of the individual components of 
energy density Pi{a) is known, the function H{a) and thus the line element in 
equation (13) is known. (Evaluating this equation at the present epoch one can 
determine the value of k; hence it is not necessary to provide this information 
separately.) Given Hq, the current value of the Hubble parameter, one can 
construct a critical density, by the definition: 

p^ = —2- = 1.88/i^ X 10-29 gm cm"^ = 2.8 x IQ^^h^MQ Mpc^^^ 

= 1.1 X lO^/i^ eV cm-^ = 1.1 x lO'^/i^ protons cm"^ (21) 

and parameterize the energy density, Pi(ao), of different components at the 
present epoch in terms of the critical density by Pi(ao) = ^iPc- [Observations 
[49,50] give h = 0.72 ±0.03 (statistical) ±0.07 (systematic)]. It is obvious from 
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equation (20) that k = corresponds to Qtot = = 1 while Qtot > 1 and 

Qtot < 1 correspond to k = ±1. When Qtot 7^ 1, equation (20), evahiated at 
the current epoch, gives {k/a^) = i7o(i^tot ~ 1); thereby fixing the value of 
{k/al); when, Q^ot = 1, it is conventional to take — 1 since its value can be 
rescaled. 



2.3 Composition of the universe 

It is important to stress that absolutely no progress in cosmology can be made 
until a relationship between p and P is provided, say, in the form of the func- 
tions Wi{a)s. This fact, in turn, brings to focus two issues which are not often 
adequately emphasized: 

(i) If we assume that the source is made of normal laboratory matter, then the 
relationship between p and P depends on our knowledge of how the equation 
of state for matter behaves at different energy scales. This information needs 
to be provided by atomic physics, nuclear physics and particle physics. Cos- 
mological models can at best be only as accurate as the input physics about 

is; any definitive assertion about the state of the universe is misplaced, if the 
knowledge about which it is based on is itself speculative or non existent 
at the relevant energy scales. At present we have laboratory results testing 
the behaviour of matter up to about 100 GeV and hence we can, in principle, 
determine the equation of state for matter up to 100 GeV. By and large, the 
equation of state for normal matter in this domain can be taken to be that of 
an ideal fiuid with p giving the energy density and P giving the pressure; the 
relation between the two is of the form P = wp with w = for non relativistic 
matter and w = (1/3) for relativistic matter and radiation. 

(ii) The situation becomes more complicated when we realize that it is en- 
tirely possible for the large scale universe to be dominated by matter whose 
presence is undetectable at laboratory scales. For example, large scale scalar 
fields dominated either by kinetic energy or nearly constant potential energy 
could exist in the universe and will not be easily detectable at laboratory 
scales. We see from (6) that such systems can have an equation of state of 
the form P = wp with w — 1 (for kinetic energy dominated scalar field) or 
w = —1 (for potential energy dominated scalar field). While the conservative 
procedure for doing cosmology would be to use only known forms of on the 
right hand side of Einstein's equations, this has the drawback of preventing 
progress in our understanding of nature, since cosmology could possibly be 
the only testing ground for the existence of forms of which are difficult to 
detect at laboratory scales. 

One of the key issues in modern cosmology has to do with the confiict in 
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Fig. 1. Cosmic inventory of energy densities. See text for description. (Figure 
adapted from [46].) 

principle between (i) and (ii) above. Suppose a model based on conventional 
equation of state, adequately tested in the laboratory, fails to account for a 
cosmological observation. Should one treat this as a failure of the cosmological 
model or as a signal from nature for the existence of a source not seen at 
laboratory scales? There is no easy answer to this question and we will focus 
on many facets of this issue in the coming sections. 

Figure 1 provides an inventory of the density contributed by different forms 
of matter in the universe. The x-axis is actually a combination Qh"- of Q 
and the Hubble parameter h since different components are measured by dif- 
ferent techniques. (Usually n = 1 or 2; numerical values are for h = 0.7.) 
The density parameter contributed today by visible, non relativistic, baryonic 
matter in the universe is about Qb ^ (0.01 — 0.2) (marked by triangles in 
the figure; different estimates are from different sources; see for a sample of 
references [51,52,53,54,55,56,57,58,59,60]). The density parameter due to radi- 
ation is about Qr pa 2 x 10~^ (marked by squares in the figure). Unfortunately, 
models for the universe with just these two constituents for the energy density 
are in violent disagreement with observations. It appears to be necessary to 
postulate the existence of: 

• Pressure-less (w = 0) non baryonic dark matter which does not couple 
with radiation and having a density of about Qdm ~ 0.3. Since it does not 
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emit light, it is called dark matter (and marked by a cross in the figure). 
Several independent techniques like cluster mass-to-light ratios [61] baryon 
densities in clusters [62,63] weak lensing of clusters [64,65] and the existence 
of massive clusters at high redshift [66] have been used to obtain a handle 
on JIdm- These observations are all consistent with Qnr = (^^dm + ^b) ~ 
Qdm~ (0.2-0.4). 

• An exotic form of matter {cosmological constant oi something similar) with 
an equation of state p a; —p (that is, w ~ —1) having a density parameter 
of about Qa 0.7 (marked by a filled circle in the figure). The evidence for 
will be discussed in section 3. 

So in addition to Hq^ at least four more free parameters are required to de- 
scribe the background universe at low energies (say, below 50 GeV). These are 
^B-i ^R, ^BM and Qa describing the fraction of the critical density contributed 
by baryonic matter, radiation (including relativistic particles like e.g, massive 
neutrinos; marked by a cross in the figure), dark matter and cosmological 
constant respectively. The first two certainly exist; the existence of last two 
is probably suggested by observations and is definitely not contradicted by 
any observations. Of these, only is well constrained and other quantities 
are plagued by both statistical and systematic errors in their measurements. 
The top two positions in the contribution to Q are from cosmological constant 
and non baryonic dark matter. It is unfortunate that we do not have labora- 
tory evidence for the existence of the first two dominant contributions to the 
energy density in the universe. (This feature alone could make most of the 
cosmological paradigm described in this review irrelevant at a future date!) 

The simplest model for the universe is based on the assumption that each of 
the sources which populate the universe has a constant wf, then equation (20) 
becomes 



where each of these species is identified by density parameter Qj and the equa- 
tion of state characterized by Wj. The most familiar form of energy densities 
are those due to pressure-less matter with Wi = (that is, non relativistic 
matter with rest mass energy density pc^ dominating over the kinetic energy 
density, pv'^/2) and radiation with Wi = (1/3). Whenever any one component 
of energy density dominates over others, P wp and it follows from the 
equation (22) (taking A; = 0, for simphcity) that 




(22) 



p (X a 



■3(l+«)). 



^2/[3(l+«;)] 



(23) 



For example, p oc a a oc t^^"^ if the source is relativistic and p oc a ^, a oc 
t^/^ if the source is non-relativistic. 
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This result shows that the past evolution of the universe is characterized by 
two important epochs (see eg. [43,44]): (i) The first is the radiation dominated 
epoch which occurs at redshifts greater than z^q ~ (^dm/^r) ~ 10^. For 
z > Zeq the energy density is dominated by hot relativistic matter and the 
universe is very well approximated as a A; = model with a{t) oc t^^^. (ii) 
The second phase occTirs for z <tC Zeq in which the universe is dominated by 
non relativistic matter and — in some cases — the cosmological constant. 
The form of a{t) in this phase depends on the relative values of ^bm and Qa- 
In the simplest case, with Odm ~ 1, J^a = 0, Ob -C I^dm the expansion is 
a power law with a{t) oc t"^/^. (When cosmological constant dominates over 
matter, a{t) grows exponentially.) 

During all the epochs, the temperature of the radiation varies as T oc a~^. 
When the temperature falls below T ki 10^ K, neutral atomic systems form 
in the universe and photons decouple from matter. In this scenario, a relic 
background of such photons with Planckian spectrum at some non-zero tem- 
perature will exist in the present day universe. The present theory is, however, 
unable to predict the value of T at t = to; it is therefore a free parameter re- 
lated Qr oc T^. 



2.4 Geometrical features of a universe with a cosmological constant 



The evolution of the universe has different characteristic features if there exists 
sources in the universe for which (1 + 3w) < 0. This is obvious from equation 
(8) which shows that if (p + 3P) = (1 -|- 3w)p becomes negative, then the 
gravitational force of such a source (with p > 0) will be repulsive. The simplest 
example of this kind of a source is the cosmological constant with w\ = — 1- 

To see the effect of a cosmological constant let us consider a universe with 
matter, radiation and a cosmological constant. Introducing a dimensionless 
time coordinate r = Hot and writing a = aoq{T) equation (20) can be cast in 
a more suggestive form describing the one dimensional motion of a particle in 
a potential 

(24) 



where 



H h ilAQ 



(25) 



This equation has the structure of the first integral for motion of a particle 
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Fig. 2. The phase portrait of the universe, with the "velocity" of the universe 
{dq/dr) plotted against the "position" g = (1 + z)~^ for different models with 
= 2.56 X W^h"'^, h = 0.5, JInr + J^a + = 1. Curves are parameterized by 
the value of JInr = 0.1,0.2,0.3,0.5,0.8,1.0 going from bottom to top as indicated. 
(Figure adapted from [46].) 

with energy £■ in a potential V{q). For models with Q = + Qnr + ^^a = 1, 
we can take £^ = so that {dq/dr) — ^Jv{q). Figure 2 is the phase portrait 
of the universe showing the velocity (dq/dr) as a function of the position q — 
(1 + z)~^ for such models. At high redshift (small q) the universe is radiation 
dominated and q is independent of the other cosmological parameters; hence 
all the curves asymptotically approach each other at the left end of the figure. 
At low redshifts, the presence of cosmological constant makes a difference and 
- in fact - the velocity q changes from being a decreasing function to an 
increasing function. In other words, the presence of a cosmological constant 
leads to an accelerating universe at low redshifts. 

For a universe with non relativistic matter and cosmological constant, the 
potential in (25) has a simple form, varying as {—a~^) for small a and (— a^) 
for large a with a maximum in between at g = gmax = (^nr/2^^a)^^'^. This 
system has been analyzed in detail in literature for both constant cosmological 
constant [67] and for a time dependent cosmological constant [68]. A wide 
variety of explicit solutions for a{t) can be provided for these equations. We 
briefly summarize a few of them. 

• If the "particle" is situated at the top of the potential, one can obtain a 
static solution with a = d = by adjusting the cosmological constant and 
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the dust energy density and taking k — 1. This solution, 

Krit = 47rG'pNR ^ \, (26) 

was the one which originally prompted Einstein to introduce the cosmolog- 
ical constant (see section 1.2). 

• The above solution is, obviously, unstable and any small deviation from the 
equilibrium position will cause a ^ or a — > oo. By fine tuning the values, 
it is possible to obtain a model for the universe which "loiters" around 
O' = ^max for a large period of time [69,70,71,24,25,26]. 

• A subset of models corresponds to those without matter and driven entirely 
by cosmological constant A. These models have k — (— 1,0, +1) and the 
corresponding expansion factors being proportional to [smh{Ht) , exp{Ht) , 
cosh(ift)] with A^ = 3H^. These line elements represent three different 
characterizations of the de Sitter spacctimc. The manifold is a four dimen- 
sional hyperboloid embedded in a fiat, five dimensional space with signature 
(H ). We shall discuss this in greater detail in section 9. 

• It is also possible to obtain solutions in which the particle starts from a — 
with an energy which is insufficient for it to overcome the potential barrier. 
These models lead to a universe which collapses back to a singularity. By 
arranging the parameters suitably, it is possible to make a{t) move away or 
towards the peak of the potential (which corresponds to the static Einstein 
universe) asymptotically [67]. 

• In the case of Qnr + ^^a = 1, the exphcit solution for a{t) is given by 

/ 3 \ 

a{t) oc ( sinh -Htj ; A; = 0; 3//^ = A (27) 

This solution smoothly interpolates between a matter dominated universe 
a{t) oc t^/^ at early stages and a cosmological constant dominated phase 
a{t) oc exp{Ht) at late stages. The transition from deceleration to acceler- 
ation occurs at ^acc — (SJIa/JInr)''^'^^ — 1, while the energy densities of the 
cosmological constant and the matter are equal at = (^^a/^^nr)^^^ — 1- 

The presence of a cosmological constant also affects other geometrical param- 
eters in the universe. Figure 3 gives the plot of dA{z) and dL{z); (note that 
angular diameter distance is not a monotonic function of z). Asymptotically, 
for large z, these have the limiting forms, 

dA{z) ^ 2{Hon^K)-'z-'; dL{z) ^ 2{Hon^j,)-h (28) 

The geometry of the spacetime also determines the proper volume of the uni- 
verse between the redshifts z and z + dz which subtends a solid angle dfl in 
the sky. If the number density of sources of a particular kind (say, galaxies. 
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Fig. 3. The left panel gives the angular diameter distance in units of cHq^ as a 
function of redshift. The right panel gives the luminosity distance in units of cH^^ 
as a function of redshift. Each curve is labelled by (J^nr, ^^a)- 




Fig. 4. The figure shows (dN/dQdz): it is assumed that n{z) = no(l + zf. The 
y-axis is in units of hoHq^. Each curve is labelled by (r^NRj^^A)- 

quasars, ...) is given by n{z), then the number count of sources per unit solid 
angle per redshift interval should vary as 



dN 
(Kldz 



n{z) 



dV n{z)alrl^{z)H-\z) 



dQdz 



(i + zy 



(29) 



Figure 4 shows (dN/d,Q.dz); it is assumed that n{z) — no{l + z)^. The y-axis 
is in units of UqHq^. 



3 Evidence for a non-zero cosmological constant 



There are several cosmological observations which suggests the existence of 
a non zero cosmological constant with different levels of reliability. Most of 
these determine either the value of Qnr or some combination of Qnr and Qa- 
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When combined with the strong evidence from the CMBR observations that 
the Qtot = 1 (see section 6) or some other independent estimate of Q^^, one 
is led to a non zero value for Q\. The most reliable ones seem to be those 
based on high redshift supernova [72,73,74] and structure formation models 
[75,76,77]. We shall now discuss some of these evidence. 

3. 1 Observational evidence for accelerating universe 

Figure 2 shows that the evolution of a universe with Vtjy ^ changes from 
a decelerating phase to an accelerating phase at late times. If H{a) can be 
observationally determined, then one can check whether the real universe had 
undergone an accelerating phase in the past. This, in turn, can be done if 
dL{z), say, can be observationally determined for a class of sources. Such a 
study of several high redshift supernova has led to the data which is shown in 
figures 5, 9. 

Bright supernova explosions are brief explosive stellar events which are broadly 
classified as two types. Type-la supernova occurs when a degenerate dwarf 
star containing CNO enters a stage of rapid nuclear burning cooking iron 
group elements (see eg., chapter 7 of [78]). These are the brightest and most 
homogeneous class of supernova with hydrogen poor spectra. An empirical 
correlation has been observed between the sharply rising light curve in the 
initial phase of the supernova and its peak luminosity so that they can serve 
as standard candles. These events typically last about a month and occurs 
approximately once in 300 years in our galaxy. (Type 11 supernova, which 
occur at the end of stellar evolution, are not useful as standard candles.) 

For any supernova, one can directly observe the apparent magnitude m [which 
is essentially the logarithm of the flux F observed] and its redshift. The abso- 
lute magnitude M of the supernova is again related to the actual luminosity 
L of the supernova in a logarithmic fashion. Hence the relation F = (L/Andl) 
can be written as 



The numerical factors arise from the astronomical conventions used in the 
definition of m and M. Very often, one will use the dimensionless combination 
{HodL{z) / c) rather than dL{z) and the above equation will change to m{z) — 
M. + 51ogiQ(i?oo?L(^)/c) with the quantity M. being related to M by 




(30) 



= M + 25 + 5 log 



( 



IMpc 



) 



= M - 5 logio h + 42.38 



(31) 
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If the absolute magnitude of a class of Type I supernova can be determined 
from the study of its light curve, then one can obtain the for these supernova 
at different redshifts. (In fact, we only need the low- 2; apparent magnitude at 
a given z which is equivalent to knowing M..) Knowing (II, one can determine 
the geometry of the universe. 

To understand this effect in a simple context, let us compare the luminosity 
distance for a matter dominated model (I^nr = 1, = 0) 

dL = [(1 + + z)^/^] , (32) 

with that for a model driven purely by a cosmological constant (JInr — 0,D,\ — 
1): 

dL^H^^z{l + z). (33) 



It is clear that at a given z, the di is larger for the cosmological constant model. 
Hence, a given object, located at a fixed redshift, will appear brighter in the 
matter dominated model compared to the cosmological constant dominated 
model. Though this sounds easy in principle, the statistical analysis turns out 
to be quite complicated. 

The supernova cosmology project (SCP) has discovered [74] 42 supernova in 
the range (0.18 — 0.83). The high- 2; supernova search team (HSST) discovered 
14 supernova in the redshift range (0.16 — 0.62) and another 34 nearby super- 
nova [73] and used two separate methods for data fitting. (They also included 
two previously published results from SCP.) Assuming + J^a = 1, the 
analysis of this data gives Qnr = 0.28 ± 0.085 (stat) ±0.05 (syst). 

Figure 5 shows the diiz) obtained from the supernova data and three theoret- 
ical curves all of which are ior k — models containing non relativistic matter 
and cosmological constant. The data used here is based on the redshift mag- 
nitude relation of 54 supernova (excluding 6 outliers from a full sample of 60) 
and that of SN 1997ff at z = 1.755; the magnitude used for SN 1997ff has been 
corrected for lensing effects [79]. The best fit curve has ^nr ~ 0.32, 0,\ 0.68. 
In this analysis, one had treated Qnr and the absolute magnitude M as free 
parameters (with -|- = 1) and has done a simple best fit for both. The 
corresponding best fit value for is = 23.92 ± 0.05. Frame (a) of figure 
6 shows the confidence interval (for 68 %, 90 % and 99 %) in the ^nr — -M. 
for the flat models. It is obvious that most of the probability is concentrated 
around the best fit value. We shall discuss frame (b) and frame (c) later on. 
(The discussion here is based on [80].) 

The confidence intervals in the — ^nr plane are shown in figure 7 for 
the full data. The confidence regions in the top left frame are obtained after 
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Fig. 5. The luminosity distance of a set of type la supernova at different red- 
sliifts and three theoretical models with J^nr + ^^A = 1- The best fit curve has 
Onr = 0.32, = 0.68. 
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Fig. 6. Confidence contours corresponding to 68 %, 90 % and 99 % based on SN 
data in the f^NR — -M plane for the flat models with I^nr + ^^A = 1- Frame (a) on 
the left uses all data while frame (b) in the middle uses low redshift data and the 
frame (c) in the right uses high redshift data. While neither the low-z or high-z data 
alone can exclude the JInr = li^^A = model, the full data excludes it to a high 
level of significance. 

marginalizing over M.. (The best fit value with la error is indicated in each 
panel and the confidence contours correspond to 68 %, 90 % and 99 %.) 
The other three frames show the corresponding result with a constant value 
for M. rather than by marginalizing over this parameter. The three frames 
correspond to the mean value and two values in the wings of Icr from the 
mean. The dashed line connecting the points (0,1) and (1,0) correspond to a 
universe with f^N^ + f^A = 1. From the figure we can conclude that: (i) The 
results do not change significantly whether we marginalize over M. or whether 
we use the best fit value. This is a direct consequence of the result in frame 
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(a) of figure (6) which shows that the probability is sharply peaked, (ii) The 
results exclude the ^nr = Ij^a = model at a high level of significance in 
spite of the uncertainty in M. 




0.5 1 1.5 2 0.5 1 1.5 2 



"NR NR 

Fig. 7. Confidence contours corresponding to 68 %, 90 % and 99 % based on SN 
data in the Onr — ^^A plane. The top left frame is obtained after marginalizing 
over Ai while the other three frames uses fixed values for A4. The values are 
chosen to be the best-fit value for Ai and two others in the wings of la limit. 
The dashed line corresponds to the flat model. The unbroken slanted line corre- 
sponds to H^dLiz = 0.63) = constant. It is clear that: (i) The data excludes the 
^NR = Ijf^A = model at a high significance level irrespective of whether we 
marginalize over M or use an accepted la range of values for ^A. (ii) The shape 
of the confidence contours are essentially determined by the value of the luminosity 
distance at z « 0.6. 

The slanted shape of the probability ellipses shows that a particular linear 
combination of ^nr and VL^ is selected out by these observations [81]. This 
feature, of course, has nothing to do with supernova and arises purely because 
the luminosity distance di depends strongly on a particular linear combination 
of and i^NR, as illustrated in figure 8. In this figure, I^nr, ^^a are treated 
as free parameters [without the = constraint] but a particular linear com- 
bination q = (0.8f2NR — 0.6fl\) is held fixed. The di is not very sensitive to 
individual values of ^nr, ^a at low redshifts when (0.8f2NR — 0.6fl\) is in the 
range of (—0.3, —0.1). Though some of the models have unacceptable parame- 
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ter values (for other reasons), supernova measurements alone cannot rule them 
out. Essentially the data at 2; < 1 is degenerate on the linear combination of 
parameters used to construct the variable q. The supernova data shows that 
most likely region is bounded by —0.3 < (0.8f2NR — 0.6Oa) ^ —0.1. In fig- 
ure 7 we have also over plotted the line corresponding to HodL{z — 0.63) = 
constant. The coincidence of this line (which roughly corresponds to rf/, at a 
redshift in the middle of the data) with the probability ellipses indicates that 
it is this quantity which essentially determines the nature of the result. 
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Fig. 8. The luminosity distance for a class of models with widely varying J1nr,^^A 
but with a constant value for q = (O.SOnr — CGJ^a) are shown in the figure. It is 
clear that when q is fixed, low redshift observations cannot distinguish between the 
different models even if J^nr and ^\ vary significantly. 

We saw earher that the presence of cosmological constant leads to an acceler- 
ating phase in the universe which — however — is not obvious from the above 
figures. To see this explicitly one needs to display the data in the a, vs a plane, 
which is done in figure 9. Direct observations of supernova is converted into 
dL{z) keeping M a free parameter. The is converted into dniz) assuming 
k — and using (17). A best fit analysis, keeping (M, Qnr) as free parame- 
ters now lead to the results shown in figure 9, which confirms the accelerating 
phase in the evolution of the universe. The curves which are over-plotted cor- 
respond to a cosmological model with I^nr + I^a = 1- The best fit curve has 
f^NR = 0.32, l^A = 0.68. 

In the presence of the cosmological constant, the universe accelerates at low 
redshifts while decelerating at high redshift. Hence, in principle, low redshift 
supernova data should indicate the evidence for acceleration. In practice, how- 
ever, it is impossible to rule out any of the cosmological models using low 
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Fig. 9. Observations of SN are converted into the 'velocity' a of the universe using 

a fitting function. The curves which are over-plotted corresponds to a cosmological 
model with ^^nr + ^^a = 1- The best fit curve has ^^nr = 0.32, VL\ = 0.68. 

redshift {z < 0.2) data as is evident from figure 9. On the other hand, high 
redshift supernova data alone cannot be used to estabhsh the existence of a 
cosmological constant. The data for [z > 0.2) in figure 9 can be moved verti- 
cally up and made consistent with the decelerating Q = 1 universe by choosing 
the absolute magnitude M suitably. It is the interplay between both the high 
redshift and low redshift supernova which leads to the result quoted above. 

This important result can be brought out more quantitatively along the fol- 
lowing lines. The data displayed in figure 9 divides the supernova into two 
natural classes: low redshift ones in the range < 2; < 0.25 (corresponding to 
the accelerating phase of the universe) and the high redshift ones in the range 
0.25 < 2; < 2 (in the decelerating phase of the universe). One can repeat all the 
statistical analysis for the full set as well as for the two individual low redshift 
and high redshift data sets. Frame (b) and (c) of figure 6 shows the confidence 
interval based on low redshift data and high redshift data separately. It is ob- 
vious that the f^NR = 1 model cannot be ruled out with cither of the two data 
sets! But, when the data sets are combined — because of the angular slant 
of the eUipses — they isolate a best fit region around Qnr ~ 0.3. This is also 
seen in figure 10 which plots the confidence intervals using just the high-z and 
low-z data separately. The right most frame in the bottom row is based on 
the low-z data alone (with marginalization over }A) and this data cannot be 
used to discriminate between cosmological models effectively. This is because 
the at low-z is only very weakly dependent on the cosmological parameters. 
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So, even though the acceleration of the universe is a low-z phenomenon, we 
cannot rcUably determine it using low-z data alone. The top left frame has 
the corresponding result with high-z data. As we stressed before, the fl^R = 1 
model cannot be excluded on the basis of the high-z data alone either. This is 
essentially because of the nature of probability contours seen earlier in frame 
(c) of figure 6. The remaining 3 frames (top right, bottom left and bottom 
middle) show the corresponding results in which fixed values of — rather 
than by marginalizing over Ai. Comparing these three figures with the cor- 
responding three frames in 7 in which all data was used, one can draw the 
following conclusions: (i) The best fit value for M. is now M. — 24.05 ± 0.38; 
the la error has now gone up by nearly eight times compared to the result 
(0.05) obtained using all data. Because of this spread, the results are sensitive 
to the value of M. that one uses, unlike the situation in which all data was 
used, (ii) Our conclusions will now depend on M.. For the mean value and 
lower end of A4, the data can exclude the Qnr — — model [see the 

two middle frames of figure 10]. But, for the high-end of allowed la range of 
A4, we cannot exclude the Onr — l,flA — model [see the bottom left frame 
of figure 10]. 



While these observations have enjoyed significant popularity, certain key points 
which undcrly these analysis need to be stressed. (For a sample of views which 
goes against the main stream, see [82,83]). 



• The basic approach uses the supernova type I light curve as a standard 
candle. While this is generally accepted, it must be remembered that we 
do not have a sound theoretical understanding of the underlying emission 
process. 

• The supernova data and fits are dominated by the region in the parameter 

space around (^nr, ^a) ~ (0.8, 1.5) which is strongly disfavoured by several 
other observations. If this disparity is due to some other unknown systematic 
effect, then this will have an effect on the estimate given above. 

• The statistical issues involved in the analysis of this data to obtain best fit 
parameters are non trivial. As an example of how the claims varied over 
time, let us note that the analysis of the first 7 high rcdshift SNc la gave a 
value for I^nr which is consistent with unity: Qnr = (0-94^g;28)- However, 
adding a single z = 0.83 supernova for which good HST data was available, 
lowered the value to Onr = (0.6 ± 0.2). More recently, the analysis of a 
larger data set of 42 high redshift SNe la gives the results quoted above. 
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Fig. 10. Confidence contours corresponding to 68 %, 90 % and 99 % based on SN 
data in the Onr — i}\ plane using either the low-z data (bottom right frame) or 
high-z data (the remaining four frames). The bottom right and the top left frames 
are obtained by marginalizing over A4 while the remaining three uses fixed values 
for Ai. The values are chosen to be the best-fit value for M and two others in the 
wings of la limit. The dashed line corresponds to the flat model. The unbroken 
slanted line corresponds to Hodiiz = 0.63) = constant. It is clear that: (i) The la 
error in top left frame (0.38) has gone up by nearly eight times compared to the 
value (0.05) obtained using all data (see figure 7) and the results are sensitive to 
the value of M. (ii) The data can exclude the I^nr = 1, = model if the mean 
or low-end value of M is used [see the two middle frames]. But, for the high-end 
of allowed la range of A^, we cannot exclude the Onr = 1, = model [see the 
bottom left frame], (iii) The low-z data [bottom right] cannot exclude any of the 
models. 

3.2 Age of the universe and cosmological constant 

From equation (24) we can also determine the current age of the universe by 
the integral 

H,to = j ^=== (34) 



Since most of the contribution to this integral comes from late times, we can 
ignore the radiation term and set VLr k, 0. When both and f^A are present 
and are arbitrary, the age of the universe is determined by the integral 
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(1 + z)^JflNK{l + zf + nA 

^ ^ (0.7 f^NR - O.Sf^A + 0.3)"°-^ 



(35) 



The integral, which cannot be expressed in terms of elementary functions, is 
well approximated by the numerical fit given in the second line. Contours of 
constant HqIq based on the (exact) integral are shown in figure 11. It is obvious 
that, for a given I^nr, the age is higher for models with fl\ ^ 0. 
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Fig. 11. Lines of constant i^o^o in the Onr — f^A plane. The eight Unes are for 
Hoto = (1.08,0.94,0.9,0.85,0.82,0.8,0.7,0.67) as shown. The diagonal line is the 
contour for models with riNR + i^A = 1- 

Observationally, there is a consensus [49,50] that h ^ 0.72 ± 0.07 and to 
13.5 ±1.5 Gyr [84]. This will give Hoto = 0.94 ± 0.14. Comparing this result 
with the fit in (35), one can immediately draw several conclusions: 

• If fljsiR > 0.1, then fl\ is non zero if Hoto > 0.9. A more reasonable assump- 
tion of Qnr > 0.3 we will require non zero Q,\ if Hoto > 0.82. 

• If we take Qnr — l,flA — and demand to > 12 Gyr (which is a conservative 
lower bound from stellar ages) will require h < 0.54. Thus a purely matter 
dominated Q = 1 universe would require low Hubble constant which is 
contradicted by most of the observations. 

• An open model with Qnr ~ 0.2, JIa = will require Hoto ~ 0.85. This still 
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requires ages on the lower side but values like h ~ 0.6, to ~ 13.5 Gyr are 
acceptable within error bars. 
• A straightforward interpretation of observations suggests maximum likeli- 
hood for Hoto = 0.94. This can be consistent with a, fl — 1 model only if 
Qnr ~ 0.3, Qa 0.7. 

If the universe is populated by dust-like matter (with w = 0) and another 
component with an equation of state parameter wx, then the age of the uni- 
verse will again be given by an integral similar to the one in equation (35) 
with replaced by ^xil + This wiU give 



Hot: 
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^Jdq 



dz 



+ z)^n^n{i + zy + nx{i + zf^'^' 
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1/2 



(36) 



The integrand varies from to (JInr + ^x)~^^^ in the range of integration for 
w < with the rapidity of variation decided by w. As a result, -f^o^o increases 
rapidly as w changes from to —3 or so and then saturates to a plateau. Even 
an absurdly negative value for w like w = —100 gives Hoto of only about 1.48. 
This shows that even if some exotic dark energy is present in the universe 
with a constant, negative w, it cannot increase the age of the universe beyond 
about Hoto ~ 1.48. 

The comments made above pertain to the current age of the universe. It is 
also possible to obtain an expression similar to (34) for the age of the universe 
at any given redshift z 

/ V f dz' , , , Hiz) 



and use it to constrain Vt\. For example, the existence of high redshift galax- 
ies with evolved stellar population, high redshift radio galaxies and age dating 
of high redshift QSOs can all be used in conjunction with this equation to 
put constrains on VL^ [85,86,87,88,89,90]. Most of these observations require 
either Vt^ 7^ or f2tot < 1 if they have to be consistent with h > 0.6. Un- 
fortunately, the interpretation of these observations at present requires fairly 
complex modeling and hence the results are not water tight. 
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3.3 Gravitational lensing and the cosmological constant 



Consider a distant source at redshift z which is lensed by some intervening 
object. The lensing is most effective if the lens is located midway between 
the source and the observer (sec, eg., page 196 of [46]). This distance will be 
(rcm/2) if the distance to the source is rem- (To be rigorous, one should be using 
angular diameter distances rather than rem for this purpose but the essential 
conclusion does not change.) To see how this result depends on cosmology, let 
us consider a source at redshift z — 2, and a lens, located optimally, in: (a) 
Q = 1 matter dominated universe, (b) a very low density matter dominated 
universe in the limit of f2 0, (c) vacuum dominated universe with fl^ — fltot- 
In case (a), dn = H~^{z) oc (1 + z)^'^^'^, so that 



For z = 2, this gives z^ = 0.608. For case (b), a oc t giving dfj oc (1 + z)^^ and 
rem(-2) oc ln(l + z). The equation to be solved is (1 + zj^) = (1 + zY^'^ which 
gives zl — 0.732 for z — 2. Finally, in the case of (c), dn is a constant giving 
rem(-s) oc z and zl — {l/2)z. Clearly, the lens redshift is larger for vacuum 
dominated universe compared to the matter dominated universe of any fl. 
When one considers a distribution of lenses, this will affect the probability 
for lensing in a manner which depends on From the observed statistics of 
lensing, one can put a bound on 

More formally, one can compute the probability for a source at redshift Zs being 
lensed in a f^A + ^nr = 1 universe (relative to the corresponding probability in 
a r^NR = 1, = model). This relative probabihty is nearly five times larger 
at Zg — 1 and about thirteen times larger for Zs — 2 in a, purely cosmological 
constant dominated universe [91,92,93,94,95,2,3]. This effect quantifies the fact 
that the physical volume associated with unit redshift interval is larger in 
models with cosmological constant and hence the probability that a light ray 
will encounter a galaxy is larger in these cases. 

Current analysis of lensing data gives somewhat differing constraints based 
on the assumptions which are made [96,97,98,99]; but typically all these con- 
straints are about fl\ < 0.7. The result [100] from Cosmic Lens All Sky Survey 
(CLASS), for example, gives Qnr = 0.311I]:?! (68%) to,\l (systematic) for a 




(38) 



The lens redshift is determined by the equation 




(39) 
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A; = universe. 



3.4 Other geometrical tests 

The existence of a maximum for dA{z) is a generic feature of cosmological 
models with Qnr > 0. For a k — 0, Qnr = 1 niodel, the maximum occurs 
at Zuiax ~ 1-25 and z^ax increases as Qa is increased. To use this as a cos- 
mological test, we require a class of objects with known transverse dimension 
and rcdshift. The most reliable quantity used so far corresponds to the phys- 
ical wavelength acoustic vibrations in the baryon-photon gas at 2; ~ 10^. 
This length scale is imprinted in the temperature anisotropies of the CMBR 
and the angular size of these anisotropies will depend on dA and hence on 
the cosmological parameters; this is discussed in section 6. In principle, one 
could also use angular sizes of galaxies, clusters of galaxies, or radio galaxies 
[101,102,103]. Unfortunately, understanding of different physical effects and 
the redshift evolution of these sources make this a difficult test in practice. 

There is another geometrical feature of the universe in which angular diameter 
distance plays an interesting role. In a closed Friedmann model with k = +1, 
there is possibihty that an observer at % = will be able to receive the light 
from the antipodal point % = tt. In a purely matter dominated universe, it is 
easy to see that the light ray from the antipodal point x = tt reaches x = 
exactly at the time of maximum expansion; therefore, in a closed, matter 
dominated universe, in the expanding phase, no observer can receive light 
from the antipodal point. The situation, however, is different in the presence 
of cosmological constant. In this case, dA{z) oc {1 + z)~^ sin /j, where 



It follows that d^ ^ when jj, ^ n. Therefore, the angular size of an object 
near the antipodal point can diverge making the object extremely bright in 
such a universe. Assuming that this phenomena does not occur up to, say 
z — 6, will imply that the redshift of the antipodal point Za{^A, ^nr) is larger 
than 6. This result can be used to constrain the cosmological parameters 
[104,105,68] though the limits obtained are not as tight as some of the other 
tests. 

Another test which can be used to obtain a handle on the geometry of the 
universe is usually called Alcock-Paczynski curvature test [106]. The basic idea 
is to use the fact that when any spherically symmetric system at high redshift 
is observed, the cosmological parameters enter differently in the characteri- 
zation of radial and transverse dimensions. Hence any system which can be 




(40) 
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approximated a priori to be intrinsically spherical will provide a way of deter- 
mining cosmological parameters. The correlation function of SDSS luminous 
red galaxies seems to be promising in terms of both depth and density for ap- 
plying this test (see for a detailed discussion, [107,108]). The main sources of 
error arises from non linear clustering and the bias of the red galaxies, either 
of which can be a source of systematic error. A variant of this method was 
proposed using observations of Lyman-a forest and compare the correlation 
function along the line of sight and transverse to the line of sight. In addi- 
tion to the modeling uncertainties, successful application of this test will also 
require about 30 close quasar pairs [109,110]. 



4 Models with evolving cosmological "constant" 

The observations which suggest the existence of non-zero cosmological con- 
stant — discussed in the last section — raises serious theoretical problems 
which we mentioned in section 1.1. These difficulties have led people to con- 
sider the possibility that the dark energy in the universe is not just a cos- 
mological constant but is of more complicated nature, evolving with time. Its 
value today can then be more naturally set by the current expansion rate 
rather than predetermined earlier on — thereby providing a solution to the 
cosmological constant problems. 

Though a host of models have been constructed based on this hope, none of 
them provides a satisfactory solution to the problems of fine-tuning. Moreover, 
all of them involve an evolving equation of state parameter wx{cl) for the 
unknown ("X") dark energy component, thereby taking away all predictive 
power from cosmology [166]. Ultimately, however, this issue needs to settled 
observationally by checking whether wx{a) is a constant [equal to —1, for the 
cosmological constant] at all epochs or whether it is indeed varying with a. 
We shall now discuss several observational and theoretical issues connected 
with this theme. 

While the complete knowledge of the [that is, the knowledge of the right 
hand side of (20)] uniquely determines H{a), the converse is not true. If we 
know only the function H(a), it is not possible to determine the nature of 
the energy density which is present in the universe. We have already seen 
that geometrical measurements can only provide, at best, the functional form 
of H{a). It follows that purely geometrical measurements of the Priedmann 
universe will never allow us to determine the material content of the universe. 

[The only exception to this rule is when we assume that each of the components 
in the universe has constant Wi. This is fairly strong assumption and, in fact, 
will allow us to determine the components of the universe from the knowledge 
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of the function H{a). To see this, we first note that the term {k/a^) in equation 
(22) can be thought of as contributed by a hypothetical species of matter with 
w = —(1/3). Hence equation (22) can be written in the form 



with a term having Wi — —(1/3) added to the sum. Let a = 3{1 + w) and 
Q{a) denote the fraction of the critical density contributed by matter with 
w = («/3) — 1. (For discrete values of Wi and ttj, the function Q{a) will be a 
sum of Dirac delta functions.) In the continuum limit, equation (41) can be 
rewritten as 

oo 

= J da n{a) e-°« (42) 

— oo 



where (a/ao) — exp(g). The function Q,{a) is assumed to have finite support 
(or decrease fast enough) for the expression on the right hand side to converge. 
If the observations determine the function H{a), then the left hand side can 
be expressed as a function of q. An inverse Laplace transform of this equation 
will then determine the form of fl{a) thereby determining the composition of 
the universe, as long as all matter can be described by an equation of state of 
the form pj = WiPi with Wi — constant for alH = 1, N.] 

More realistically one is interested in models which has a comphcated form of 
wxia) for which the above analysis is not applicable. Let us divide the source 
energy density into two components: Pk{(i), which is known from independent 
observations and a component px(a) which is not known. From (20), it follows 
that 

^p.(a) = i7»(l-Q(a)); W - (43) 

Taking a derivative of In px(a) and using (19), it is easy to obtain the relation 

^x(a) = -^-^ln[(l - Q{a))H'{a)a'] (44) 
3 a ma 

If geometrical observations of the universe give us H{a) and other observa- 
tions give us Pk{a) then one can determine Q and thus Wxia)- While this is 
possible, in principle the uncertainties in measuring both H and Q makes this 
a nearly impossible route to follow in practice. In particular, one is interested 
in knowing whether w evolves with time or a constant and this turns out to 
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be a very difficult task observationally. We shall now briefly discuss some of 
the issues. 

4.1 Parametrized equation of state and cosmological observations 

One simple, phenomenological, procedure for comparing observations with 
theory is to parameterize the function w{a) in some suitable form and deter- 
mine a finite set of parameters in this function using the observations. Theo- 
retical models can then be reduced to a finite set of parameters which can be 
determined by this procedure. To illustrate this approach, and the difficulties 
in determining the equation of state of dark energy from the observations, we 
shall assume that w{a) is given by the simple form: w{a) = Wq + Wi{l — a); 
in the k = model (which we shall assume for simplicity), Wq measures the 
current value of the parameter and —wi gives its rate of change at the present 
epoch. In addition to simplicity, this parameterization has the advantage of 
giving finite w in the entire range < a < 1. 
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Fig. 12. Confidence interval contours in the wq — wi plane arising from the full 
supernova data, for flat models with Onr + ^^A = 1- The three frames are 
for r^NR = (0.2,0.3,0.4). The data cannot rule out cosmological constant with 
wq = —l,wi = 0. The slanted line again corresponds to HQdi{z = 0.63) = con- 
stant and shows that the shape of the probability ellipses arises essentially from this 
feature. 

Figure 12 shows confidence interval contours in the Wq — Wi plane arising from 
the full supernova data, obtained by assuming that fl^^ -|- f^A = 1. The three 
frames are for Qtssr = (0.2,0.3,0.4). The following features are obvious from 
the figure: (i) The cosmological constant corresponding to Wq = —l,Wi = 
is a viable candidate and cannot be excluded. (In fact, different analysis of 
many observational results lead to this conclusion consistently; in other words, 
at present there is no observational motivation to assume Wi 7^ 0.) (ii) The 
result is sensitive to the value of I^nr which is assumed. This is understandable 
from equation (44) which shows that wx{ci) depends on both Q oc f2NR and 
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H{a). (We shall discuss this dependence of the results on J^nr in greater 
detail below), (iii) Note that the axes are not in equal units in figure 12. The 
observations can determine Wq with far greater accuracy than Wi. (iv) The 
slanted fine again corresponds to HodL{z = 0.63) = constant and shows that 
the shape of the probabihty eUipses arises essentially from this feature. 

In summary, the current data definitely supports a negative pressure com- 
ponent with Wo < —(1/3) but is completely consistent with wi = 0. If this 
is the case, then the cosmological constant is the simplest candidate for this 
negative pressure component and there is very little observational motiva- 
tion to study other models with varying w{a). On the other hand, the cos- 
mological constant has well known theoretical problems which could possi- 
bly be alleviated in more sophisticated models with varying w{a). With this 
motivation, there has been extensive amount of work in the last few years 
investigating whether improvement in the observational scenario will allow 
us to determine whether Wi is non zero or not. (For a sample of references, 
see [111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126].) In the 
context of supernova based determination of di, it is possible to analyze the 
situation along the following lines [80] . 

Since the supernova observations essentially measure dL^a), accuracy in the 
determination of Wq and Wi from (both the present and planned future [127]) 
supernova observations will crucially depend on how sensitive di is to the 
changes in wq and Wi. A good measure of the sensitivity is provided by the 
two parameters 

A{z, Wo, wi) = -f— \n{dL{z, wq, wi)Ho); 
dwo 

B{z, Wq, Wi) = l-a{dL{z, Wq, Wi)Ho) (45) 
dwi 

Since dL{z, Wq, Wi) can be obtained from theory, the parameters A and B can 
be computed form theory in a straight forward manner. At any given redshift 
z, we can plot contours of constant A and B in the wq — wi plane. Figure 
(13) shows the result of such an analysis [80]. The two frames on the left are 
at 2; = 1 and the two frames on the right are at 2; = 3. The top frames give 
contours of constant A and bottom frame give contours of constant B. From 
the definition in the equation (45) it is clear that A and B can be interpreted 
as the fractional change in di for unit change in wo,Wi. For example, along 
the line marked A — 0.2 (in the top left frame) di will change by 20 per cent 
for unit change in wq. It is clear from the two top frames that for most of 
the interesting region in the Wq — Wi plane, changing Wq by unity changes 
di hj about 10 per cent or more. Comparison of z = 1 and z = 3 (the 
two top frames) shows that the sensitivity is higher at high redshift, as to 
be expected. The shaded band across the picture corresponds to the region 
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in which —l<w{a)<0 which is of primary interest in constraining dark 
energy with negative pressure. One concludes that determining wq from 
fairly accurately will not be too daunting a task. 
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Fig. 13. Sensitivity of di to the parameters wo,wi. The curves correspond to con- 
stant values for the percentage of change in (IlHq for unit change in wq (top frames), 
and wi (bottom frames). Comparison of the top and bottom frames shows that (IlHq 
varies by few tens of percent when wq is varied but changes by much lesser amount 
whenuji is varied. 

The situation, however, is quite different as regards wi as illustrated in the 
bottom two frames. For the same region of the wq — wi plane, di changes only 
by a few percent when wi changes by unity. That is, di is much less sensitive 
to Wi than to Wq. It is going to be significantly more difficult to determine a 
value for wi from observations of d^ in the near future. Comparison of 2; = 1 
and z = 3 again shows that the sensitivity is somewhat better at high redshifts 
but only marginally. 

The situation is made worse by the fact that d^ also depends on the parameter 
^NR. If varying fl^n mimics the variation of wi or wo, then, one also needs to 
determine the sensitivity of di to Q^r. Figure 14 shows contours of constant 
Hod^ in the — wq and ^nr — Wi planes at two redshifts z = 1 and z = 3. 
The two top frames shows that if one varies the value of ^nr in the allowed 
range of, say, (0.2, 0.4) one can move along the curve of constant d^ and induce 
fairly large variation in wi. In other words, large changes in wi can be easily 
compensated by small changes in ^nr while maintaining the same value for di 
at a given redshift. This shows that the uncertainty in Qnr introduces further 
difficulties in determining wi accurately from measurements of d^. The two 
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Fig. 14. Contours of constant HqcIl in the f^NR — wq and Onr — wi planes at two 
redshifts z = 1 and z = 3. The two top frames shows that a smah variation of 
JInr in the allowed range of, say, (0.2, 0.4)corresponds to fairly large variation in 
wi along the curve of constant d-L- 

lower frames show that the situation is better as regards Wq. The curves are 
much less steep and hence varying J^nr does not induce large variations in wq- 
We are once again led to the conclusion that unambiguous determination of 
Wi from data will be quite difficult. This is somewhat disturbing since wi ^ 
is a clear indication of a dark energy component which is evolving. It appears 
that observations may not be of great help in ruling out cosmological constant 
as the major dark energy component. (The results given above are based on 
[80]; also see [128] and references cited therein.) 

4-2 Theoretical models with time dependent dark energy: cosmic degeneracy 

The approach in the last section was purely phenomenological and one might 
like to construct some physical model which leads to varying w{a). It turns 
out that this is fairly easy, and — in fact — - it is possible to construct models 
which will accommodate virtually any form of evolution. We shall now discuss 
some examples. 

A simple form of the source with variable w are scalar fields with Lagrangians 
of different forms, of which we will discuss two possibilities: 

L^uin = lda<l)dy - V{<j)y, Ltach = -Vmi - da(l)d^(t>f' (46) 
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Both these Lagrangians involve one arbitrary function V{(f)). The first one, 
-Squill, which is a natural generalisation of the Lagrangian for a non-relativistic 
particle, L = (l/2)g^ — V{q), is usually called quintessence (for a sample of 
models, see [129,130,131,132,133,134,135,136,137,138,139]). When it acts as a 
source in Friedman universe, it is characterized by a time dependent w{t) with 



The structure of the second scalar field can be understood by a simple anal- 
ogy from special relativity. A relativistic particle with (one dimensional) po- 
sition q{t) and mass m is described by the Lagrangian L = —m\J\ — . 
It has the energy E = mj \/\ — q^ and momentum p = mqj \/\ — q^ which 
are related by E'^ — + m^. As is well known, this allows the possibihty 
of having massless particles with finite energy for which E"^ = p^. This is 
achieved by taking the limit of m ^ and g — > 1, while keeping the ratio in 
E = m/ yjl — q^ finite. The momentum acquires a life of its own, unconnected 
with the velocity g, and the energy is expressed in terms of the momentum 
(rather than in terms of q) in the Hamiltonian formulation. We can now con- 
struct a field theory by upgrading q{t) to a field 0. Relativistic invariance 
now requires </> to depend on both space and time [(/> = 0(t, x)] and q^ to be 
replaced by di(f)d^(f). It is also possible now to treat the mass parameter m 
as a function of 0, say, V{(j)) thereby obtaining a field theoretic Lagrangian 
L — —V{(f))\/1 — d^(pdi(p. The Hamiltonian structure of this theory is alge- 
braically very similar to the special relativistic example we started with. In 
particular, the theory allows solutions in which V ^ 0, di(j)d^(j) ~^ 1 simul- 
taneously, keeping the energy (density) finite. Such solutions will have finite 
momentum density (analogous to a massless particle with finite momentum 
p) and energy density. Since the solutions can now depend on both space 
and time (unlike the special relativistic example in which q depended only on 
time), the momentum density can be an arbitrary function of the spatial coor- 
dinate. This provides a rich gamut of possibilities in the context of cosmology. 

[140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164 
This form of scalar field arises in string theories [167] and — for technical rea- 
sons — is called a tachyonic scalar field. (The structure of this Lagrangian is 
similar to those analyzed in a wide class of models called K-essence; see for 
example, [160]. We will not discuss K-essence models in this review.) 

The stress tensor for the tachyonic scalar field can be written in a perfect fiuid 
form 



1 - i2V/ct>') 
1 + (2^/02) 



(47) 



Ti = {p + p)u'uk - pSi 



(48) 
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with 



"' = 7^^ "=7^?^ P^-VWyrr^ (49) 

The remarkable feature of this stress tensor is that it could be considered as 

the sum of a pressure less dust component and a cosmological constant [165] 
To show this explicitly, wc break up the density p and the pressure p and write 
them in a more suggestive form as 

P = Pa + Pdm; p^Pv+ Pdm (50) 



where 



Pdm = ^i^l^J'^g^f ^ °' ^(<^)Vl - Pv = -Pa(51) 



This means that the stress tensor can be thought of as made up of two com- 
ponents - one behaving like a pressure-less fluid, while the other having a 
negative pressure. In the cosmological context, the tachyonic field is described 
by: 

p,{t) = v[i - 02]-'/'; Pi = -v[i - <pV; Wi = ^'-i (52) 



When <j) is small (compared to V in the case of quintessence or compared to 
unity in the case of tachyonic field), both these sources have ty — > — 1 and 
mimic a cosmological constant. When (p ^ V, the quintessence has w fa 1 
leading to pg (x (1 + z)^; the tachyonic field, on the other hand, has w ~ for 
(p ^ 1 and behaves like non-relativistic matter. In both the cases, — 1 < w < 1, 
though it is possible to construct more complicated scalar field Lagrangians 
with even w < —1. (See for example, [168]; for some other alternatives to 
scalar field models, see for example, [169].) 

Since the quintessence field (or the tachyonic field) has an undetermined free 
function V^(0), it is possible to choose this function in order to produce a given 
H{a). To see this exphcitly, let us assume that the universe has two forms of 
energy density with p{a) = pknown(a) + P</)(a) where pkiiown(a) arises from any 
known forms of source (matter, radiation, ...) and p,/,(a) is due to a scalar field. 
When w{a) is given, one can determine the ^(0) using either (47) or (52). For 
quintessence, (47) along with (43) gives 



02(a) =p(l + «;) = ^1^(1 -g)(l + «;); 
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2V{a)=p{l-w) = ^^{l-Q){l-w) 
For tachyonic scalar field, (52) along with (43) gives 



(/.^(a) = (l + w); V{a)^p{-wy/^^ 



SttG 



(i-g)(-^. 



,1/2 



(53) 



(54) 



Given Q{a), w{a) these equations implicitly determine V"(0). Wc have already 
seen that, for any cosmological evolution specified by the functions H{a) and 
Pfc(a), one can determine w{a); see equation (44). Combining (44) with either 
(53) or (54), one can completely solve the problem. 

Let us first consider quintessence. Here, using (44) to express w in terms of H 
and Q, the potential is given implicitly by the form [170,166] 



V{a) 
0(a) = 



167rG 



H{l-Q) 



6H + 2aH' - 



aHQ' 



- 1 ■ 


1/2 


r da 




J 


a 



aQ' - (1 - g) 



1-Q 

din a 



1/2 



(55) 
(56) 



where Q{a) = [87rGpknown(o)/3i?^(a)]. We shall now discuss some examples of 
this result: 

• Consider a universe in which observations suggest that iJ^(a) = H'^ar^ . 
Such a universe could be populated by non relativistic matter with density 
parameter Onr = = 1. On the other hand, such a universe could be popu- 
lated entirely by a scalar field with a potential T^(0) = Vq exp[— (167rG/3)^/^0] 
One can also have a linear combination of non relativistic matter and scalar 
field with the potential having a generic form ViS) = Aexp[— 50]. 

• Power law expansion of the universe can be generated by a quintessence 
model with V{(j)) = 0~". In this case, the energy density of the scalar field 
varies as oc ^-2a/(2+Q;). jf ^j^g background density pbg varies as pbg oc 

the ratio of the two energy densities changes as (P(/,/pbg = t'^^^'^^"'^)- 
Obviously, the scalar field density can dominate over the background at 
late times for a > 0. 

• A different class of models arise if the potential is taken to be exponential 
with, say, V{(t)) oc exp(— A0/Mpi). When A; = 0, both p^ and pbg scale in the 
same manner leading to 

P4> ^ 3(1 + Wbg) ^^^^ 
Pbg + P<t> 

where Wbg refers to the background parameter value. In this case, the dark 
energy density is said to "track" the background energy density. While this 



42 



could be a model for dark matter, there are strong constraints on the total 
energy density of the universe at the epoch of nucleosynthesis. This requires 
fl(f, < 0.2 requiring dark energy to be sub dominant at all epochs. 
Many other forms of H{a) can be reproduced by a combination of non- 
relativistic matter and a suitable form of scalar field with a potential V{(l)). 
As a final example [68], suppose H'^{a) = i7Q[nNRa~^ + (1 ^ ^^nk This 
can arise, if the universe is populated with non-relativistic matter with 
density parameter fi^R and a scalar field with the potential, determined 
using equations (55), (56). We get 



V{(f>) = Hsinh'"/^"-^) [«(</) - V')] 



(58) 



where 



(6 - n)H^ 



IGttG 
and ■0 is a constant. 



"nr 



NR 



1 

n — 3 



a = {3-n){2nG/n 



,1/2 



(59) 



Similar results exists for the tachyonic scalar field as well [166]. For example, 
given any H{t), one can construct a tachyonic potential V{(f)) so that the scalar 
field is the source for the cosmology. The equations determining V{(j)) are now 
given by: 



... r da f aQ' 



Q) 




(60) 



V 



MP 
8^ 



2aH' 



aQ' 



1/2 



3 H 3(1 -Q) 



(61) 



Equations (60) and (61) completely solve the problem. Given any H{t), these 
equations determine V{t) and 0(t) and thus the potential 1^(0). 

As an example, consider a universe with power law expansion a = t". If it 
is populated only by a tachyonic scalar field, then Q — 0; further, (aH'/H) 
in equation (60) is a constant making a constant. The complete solution is 
given by 

m-Q t^*, vw = ^(i-A) ^ (62) 



where n > (2/3). Combining the two, we find the potential to be 
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For such a potential, it is possible to have arbitrarily rapid expansion with 
large n. (For the cosmological model, based on this potential, see [159].) 

A wide variety of phenomenological models with time dependent cosmologi- 
cal constant have been considered in the literature. They involve power law 
decay of cosmological constant like A oc [171,172,173,174,175,176,68] or 
A oc a-", [177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192], 
exponential decay A oc exp(— aa) [193] and more complicated models (for a 
summary, see [68]). Virtually all these models can be reverse engineered and 
mapped to a scalar field model with a suitable 1^(0). Unfortunately, all these 
models lack predictive power or clear particle physics motivation. 

This discussion also illustrates that even when w{a) is known, it is not possible 
to proceed further and determine the nature of the source. The explicit exam- 
ples given above shows that there are at least two different forms of scalar field 
Lagrangians (corresponding to the quintessence or the tachyonic field) which 
could lead to the same w{a). A theoretical physicist, who would like to know 
which of these two scalar fields exist in the universe, may have to be content 
with knowing w{a). The accuracy of the determination of w{a) depends on 
the prior assumptions made in determining Q, as well as on the observational 
accuracy with which the quantities H{a) can be measured. Direct observations 
usually give the luminosity distance (i^, or angular diameter distance d^- To 
obtain H{a) from either of these, one needs to calculate a derivative [see, for 
example, (17)] which further limits the accuracy significantly. As we saw in 
the last section, this is not easy. 



5 Structure formation in the universe 

The conventional paradigm for the formation of structures in the universe is 
based on the growth of small perturbations due to gravitational instabilities. 
In this picture, some mechanism is invoked to generate small perturbations in 
the energy density in the very early phase of the universe. These perturbations 
grow due to gravitational instability and eventually form the different struc- 
tures which we see today. Such a scenario can be constrained most severely by 
observations of cosmic microwave background radiation (CMBR) at 2; ^ 10^. 
Since the perturbations in CMBR are observed to be small (10^'"^ — 10^^ de- 
pending on angular scales), it follows that the energy density perturbations 
were small compared to unity at the redshift of z 1000. 

The central quantity one uses to describe the growth of structures during 
< z < 10^ is the density contrast defined as (5(t, x) = [p(t, x) — Phg{t)]/ phgit) 
which characterizes the fractional change in the energy density compared to 
the background. (Here Pbg(^) is the mean background density of the smooth 
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universe.) Since one is often interested in the statistical behaviour of structures 
in the universe, it is conventional to assume that 6 and other related quantities 
are elements of an ensemble. Many popular models of structure formation 
suggest that the initial density perturbations in the early universe can be 
represented as a Gaussian random variable with zero mean (that is, < 5 >— 0) 
and a given initial power spectrum. The latter quantity is defined through the 
relation P{t,k) =< > where 5k is the Fourier transform of 5(t,x) 

and < ... > indicates averaging over the ensemble. It is also conventional to 
define the two-point correlation function C,{t,x) as the Fourier transform of 
P(i, k) over k. Though gravitational clustering will make the density contrast 
non Gaussian at late times, the power spectrum and the correlation function 
continue to be of primary importance in the study of structure formation. 

When the density contrast is small, its evolution can be studied by linear 
perturbation theory and each of the spatial Fourier modes ^^(t) will grow in- 
dependently. It follows that d(t, x) will have the form d(t, x) = D(t)/(x) in the 
linear regime where D{t) is the growth factor and /(x) depends on the initial 
configuration. When 5 ~ 1, linear perturbation theory breaks down and one 
needs to either use some analytical approximation or numerical simulations to 
study the non linear growth. A simple but effective approximation is based on 
spherical symmetry in which one studies the dynamics of a spherical region in 
the universe which has a constant over-density compared to the background. 
As the universe expands, the over-dense region will expand more slowly com- 
pared to the background, will reach a maximum radius, contract and virialize 
to form a bound nonlinear system. If the proper coordinates of the particles in 
a background Friedmann universe is given by r = a(t)x we can take the proper 
coordinates of the particles in the over-dense region to be r = i?(t)x where 
R{t) is the expansion rate of the over-dense region. The relative acceleration 
of two geodesies in the over-dense region will be g = Rx = {R/R)r. Using (8) 
and V • r = 3, we get 

" AttG, GM AnG, 

R = ^(P + 3P)i? = -—^ —{p + 3P)nondusti? (64) 

where the subscript 'non-dust' refers to all components of matter other than 
the one with equation of state P = 0; the dust component is taken into account 
by the first term on the right hand side with M = (47r/3)pNR-R^- The density 
contrast is related to R by {1 + S) = (p/pbg) = (a/P)^. Given the equation 
(64) satisfied by R and (20), it is easy to determine the equation satisfied by 
the density contrast. We get (see p. 404 of [9]): 

S + 47rGp,(l + S)S + (65) 

a 3(1-1-0) 

This is a fully nonlinear equation satisfied by the density contrast in a spher- 
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ically symmetric over-dense region in the universe. 
5.1 Linear evolution of perturbations 

When the perturbations are small, one can ignore the second term in the right 
hand side of (65) and replace (1 + 5) by unity in the first term on the right 
hand side. The resulting equation is vahd in the hnear regime and hence will be 
satisfied by each of the Fourier modes 5]^{t) obtained by Fourier transforming 
(5(t,x) with respect to x. Taking x) = D{t)f{'x), the D{t) satisfies the 
equation 

D + 2-D = AnGpbD (66) 
a 

The power spectra P{k, t) =< |5k(^)P > at two different redshifts in the linear 
regime are related by 

P{k, Zf) = T'{k, Zf, Zi, hg)P{k, Zi) (67) 

where T (called transfer function) depends only on the parameters of the 
background universe (denoted by 'bg') but not on the initial power spectrum 
and can be computed by solving (66). It is now clear that the only new input 
which structure formation scenarios require is the specification of the initial 
perturbation at all relevant scales, which requires one arbitrary function of 
the wavenumber k — 27r/A. 

Let us first consider the transfer function. The rate of growth of small pertur- 
bations is essentially decided by two factors: (i) The relative magnitudes of 
the proper wavelength of perturbation Aprop(t) oc a{t) and the Hubble radius 
(i/j(t) = H^^(t) = (a/a)^^ and (ii) whether the universe is radiation domi- 
nated or matter dominated. At sufficiently early epochs, the universe will be 
radiation dominated and dnit) oc t will be smaller than the proper wavelength 
Aprop(^) t^l"^. The density contrast of such modes, which are bigger than the 
Hubble radius, will grow as until Aprop = dnit). [When this occurs, the 
perturbation at a given wavelength is said to enter the Hubble radius. One 
can use (66) with the right hand side replaced by 47r(l +'u;)(l + 3w)Gp in this 
case; this leads to oc t oc a^.] When Ap^p < dn and the universe is radia- 
tion dominated, the perturbation does not grow significantly and increases at 
best only logarithmically [194]. Later on, when the universe becomes matter 
dominated for t > t^q, the perturbations again begin to grow. It follows from 
this result that modes with wavelengths greater than dcq = c^i/(^cq) — which 
enter the Hubble radius only in the matter dominated epoch — continue to 
grow at all times while modes with wavelengths smaller than deq suffer lack 
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of growth (in comparison with longer wavelength modes) during the period 
Center < ^ < ^cq- This fact leads to a distortion of the shape of the primordial 
spectrum by suppressing the growth of small wavelength modes in comparison 
with longer ones. Very roughly, the shape of T'^{k) can be characterized by the 
behaviour T^{k) oc for k > keq and fa 1 for A; < keq. The wave number 
keq corresponds to the length scale 

4q = dnizeq) = {27r / k^q) ^ 13(1^/^2) "^Mpc (68) 



(eg., [44], p. 75). The spectrum at wavelengths A » deq is undistorted by the 

evolution since is essentially unity at these scales. Further evolution can 
eventually lead to nonlinear structures seen today in the universe. 

At late times, we can ignore the effect of radiation in solving (66). The linear 
perturbation equation (66) has an exact solution (in terms of hyper-geometric 
functions) for cosmological models with non-relativistic matter and dark en- 
ergy with a constant w. It is given by 



D(a) 



2F1 



1 w-1 ^ 5 

2w ^ 6w^ 



■3w 



1-n 



NR 



NR 



(69) 



[This result can be obtained by direct algebra. If the independent variable 
in equation (66) is changed from t to a"^'^ and the dependent variable is 
changed from D to (D/a), the resulting equation has the standard form of 
hypergeometric equation for a universe with dark energy and non-relativistic 
matter as source.] Figure 15 shows the growth factor for different values of w 
including the one for cosmological constant (corresponding to w — —1) and 
an open model (with w — —1/3.) 

For small values of a, D ^ a which is an exact result for = 0, Qnr = 1 
model. The growth rate slows down in the cosmological constant dominated 
phase (in models with fl^R + flA — ^ with w — —1) or in the curvature 
dominated phase (open models with Qnr < 1 corresponding to w — —1/3). 
Between the two cases, there is less growth in open models compared to models 
with cosmological constant. 

It is possible to rewrite equation (65) in a different form to find an approximate 
solution for even variable w{a). Converting the time derivatives into derivatives 
with respect to a (denoted by a prime) and using the Friedmann equations, 
we can write (65) as 
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Fig. 15. The growth factor for different values of w including the one for cosmological 
constant (corresponding to u; = — 1) and an open model (with w = —1/3). 

In a universe populated by only non relativistic matter and dark energy char- 
acterized by an equation of state function w{a), this equation can be recast in 
a different manner by introducing a time dependent Q [as in equation (43)] by 
the relation Q(t) = {8ttG /3)[pj^Bit) / H^(t)] so that {dQ/d\na) = 3u;Q(l-Q). 
Then equation (65) becomes in terms of the variable / = {dln6/dlna) 



3wQ(l-Q)^+f + J 



= |«(i + *) + i 



1 



(71) 



Unfortunately this equation is not closed in terms of / and Q since it also 
involves 6 = exp[/ {da/a)f]. But in the linear regime, we can ignore the second 
term on the right hand side and replace (i+S) by unity in the first term thereby 
getting a closed equation: 



3wQ{l-Q)^ + f + f 



-w{l-Q) 



2^ 



(72) 



This equation has approximate power law solutions [190] of the form f — Q"^ 
when \dw/dQ\ <^ 1/(1 — Q). Substituting this ansatz, we get 
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[Note that Q{t) — > 1 at high redshifts, which is anyway the domain of va- 
hdity of the hnear perturbation theory]. This rcsuh shows that n is weakly 
dependent on Onr; further, n ~ (4/7) for open Friedmann model with non 
relativistic matter and n (6/11) ^ 0.6 in a, k — model with cosmological 
constant. 

Let us next consider the initial power spectrum P{k, Zi) in (67). The following 
points need to be emphasized regarding the initial fluctuation spectrum. 

(1) It can be proved ihaX known local physical phenomena, arising from laws 
tested in the laboratory in a medium with (P/ p) > 0, are incapable producing 
the initial fluctuations of required magnitude and spectrum (eg., [9], p 458). 
The initial fluctuations, therefore, must be treated as arising from physics 
untested at the moment. 

(2) Contrary to claims sometimes made in the literature, inflationary models 
are not capable of uniquely predicting the initial fluctuations. It is possible 
to come up with viable inflationary potentials ([197], chapter 3) which are 
capable of producing any reasonable initial fluctuation. 

A prediction of the initial fluctuation spectrum was indeed made by Harrison [198] 
and Zeldovich [199], who were years ahead of their times. They predicted — 
based on very general arguments of scale invariance — that the initial fluc- 
tuations will have a power spectrum P = Ak^ with n = 1. Considering the 
simplicity and importance of this result, we shall briefly recall the arguments 
leading to the choice of n = 1. 

If the power spectrum is P oc A;" at some early epoch, then the power per 
logarithmic band of wave numbers is oc k^P{k) oc k^"'~^^\ Further, when 
the wavelength of the mode is larger than the Hubble radius, dnit) = {d/a)~^, 
during the radiation dominated phase, the perturbation grows as making 
oc a^A;^"'^^^ We need to determine how A scales with k when the mode 
enters the Hubble radius dnit). The epoch ttgnter which this occurs is de- 
termined by the relation 2Traenter/k = dn- Using oc i oc in the radiation 
dominated phase, we get Qenter oc k~^ so that 

A^{k, a^nter) OC al^teM-^-^^^ ^ ^^""'^ (74) 

It follows that the amplitude of fluctuations is independent of scale k at the 
time of entering the Hubble radius, only if n = 1. This is the essence of 
Harrison-Zeldovich and which is independent of the inflationary paradigm. 
It follows that verification of n = 1 by any observation is not a verifica- 
tion of inflation. At best it verifies a far deeper principle of scale invariance. 
We also note that the power spectrum of gravitational potential P^ scales as 
P^ (X. P/k^ /c^""^). Hence the fluctuation in the gravitational potential (per 
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decade in k) oc k^P,f, is proportional to oc k^'^~^\ This fluctuation in 
the gravitational potential is also independent of k for n = 1 clearly showing 
the special nature of this choice. [It is not possible to take n strictly equal 
to unity without specifying a detailed model; the reason has to do with the 
fact that scale invariance is always broken at some level and this will lead 
to a small difference between n and unity]. Given the above description, the 
basic model of cosmology is based on seven parameters. Of these 5 param- 
eters (i^o, i^B, ^^DM, ^A, ^r) determine the background universe and the two 
parameters (^4, n) specify the initial fluctuation spectrum. 

It is possible to provide simple analytic fitting functions for the transfer func- 
tion, incorporating all the above effects. For models with a cosmological con- 
stant, the transfer function is well fitted by [195] 

Tlip) = + 3.89P + (16.1p)^ + i5A6pf + (6.71p)r ^^^(75) 



where p = k/{Th Mpc""*^) and V = f2NR/iexp[— QB(l + \/2/i/f^NR)] is called the 
'shape factor'. The presence of dark energy, with a constant w, will also affect 
the transfer function and hence the final power spectrum. An approximate 
fitting formula can be given along the following lines [196]. Let the power 
spectrum be written in the form 

P{k,a)^AQk-T^{k) (76) 

\9q,oJ 



where Aq is a normalization, Tq is the modified transfer function and gq = 
(D/a) is the ratio between linear growth factor in the presence of dark energy 
compared to that in O, — 1 model. Writing Tq as the product Tq\Ta where 
Ta is given by (75), numerical work shows that 

TQA{k, a = = — , q = — — , 77 



where k is in Mpc , and a is a scale-independent but time-dependent coeffi- 
cient well approximated by a = {—wY with 



s = (0.012 - 0.036 w - 0.017/w)[l - VL^K{a)] 
+ (0.098 + 0.029 w - 0.085/'u;) InQNR(a) (78) 

where the matter density parameter is r2NR(a) = f2NR/[fiNR+ (1~^nr) a~^"'] . 
Similarly, the relative growth factor can be expressed in the form qq^ = 
{9q/9a) = {-wY with 
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t = -(0.255 + 0.305 w + 0.0027/m;)[1 - (]NR(a)] 
- (0.366 + 0.266 w - 0.07 /w) In n^Ria) (79) 

Finally the amplitude Aq can be expressed in the form Aq — 6^(0/ HqY^^^ / (47r) , 
where 

= 2 X 10-^ao ' (^^nr)'=^+^^'"''^^ exp [c^{n - 1) + c^{n - if] (80) 



and Q!o = ot{a = 1) of equation (78), and 



ci = -0.789|w|°-°^^^-°-='"^°l"l, C2 = -0.118 - 0.0727^/;, 

C3 = -1.037, C4 = -0.138 (81) 

This fit is valid for -1 < w < -0.2. 



5. 2 Nonlinear growth of perturbations 



In a purely matter dominated universe, equation (64) reduces to R = —GM/R^. 
Solving this equation one can obtain the non linear density contrast 5 as a 
function of the redshift z: 

n + - (^S" '^-(^ + ^-) - (^^ ^0 (..^ 

>~\2,) (^-sin^^)2/3 - bAsj (^-sin^)2/3' ^^^^ 



9 - sin 



2 (1 - cos ef 



- 1. (83) 



Here, (5j > is the initial density contrast at the redshift Zi and 6q is the 
density contrast at present if the initial density contrast was evolved by linear 
approximation. In general, the linear density contrast ^z, is given by 



When e = {2n/3), 6l = 0.568 and 6 = 1.01 ~ 1. If we interpret 5 = 1 as the 
transition point to nonlinearity, then such a transition occurs at 6' = (27r/3), 
Sl — 0.57. Prom (82), we see that this occurs at the redshift (1 + z^i) — 
((5o/0.57). The spherical region reaches the maximum radius of expansion at 
6 = n. This corresponds to a density contrast of 6m ~ 4.6 which is definitely 
in the nonlinear regime. The linear evolution gives 6l = 1.063 at = tt. 
After the spherical over dense region turns around it will continue to contract. 
Equation (83) suggests that at ^ = 27r all the mass will collapse to a point. 
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A more detailed analysis of the spherical model [200], however, shows that 
the virialized systems formed at any given time have a mean density which is 
typically 200 times the background density of the universe at that time in a 
^NR = 1- This occurs at a redshift of about (1+^coii) = (5o/l-686). The density 
of the virialized structure will be approximately pcoU — 170po(l + -^cou)^ where 
po is the present cosmological density. The evolution is described schematically 
in figure 16. 



9 =(2jt/3) 9=Jt 9= 2jt 




nl ^max ^coll 
Expansion Factor 



Fig. 16. Evolution of an over dense region in spherical top-hat approximation. 

In the presence of dark energy, one cannot ignore the second term in equation 
(64). In the case of a cosmological constant, w = —1 and p = constant and 
this extra term is independent of time. This allows one to obtain the first 
integral to the equation (64) and reduce the problem to quadrature (see, for 
example [201,202,203]). For a more general case of constant w with w 7^ — 1, 
the factor (p + 3P) = p(l + 3w) will be time dependent because p will be time 
dependent even for a constant w if -u; 7^ —1. In this case, one cannot obtain an 
energy integral for the equation (64) and the dynamics has to be determined 
by actual numerical integration. Such an analysis leads to the following results 
[190], [204], [205]: 

(i) In the case of matter dominated universe, it was found that the linear 
theory critical threshold for collapse, 6c, was about 1.69. This changes very 
little in the presence of dark energy and an accurate fitting function is given 

by 
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= — [1 + a logio ^nr] , 

a = 0.353w^ + 1.044^3 + 1.128^^ + 0.555w + 0.131. 




(85) 



(ii) The over density of a virialized structure as a function of tiie redsiiift of 
virialization, however, depends more sensitively on the dark energy compo- 
nent. For —l<w< —0.3, this can be fitted by the function 



a = 0.399- 1.309(|w|°-^^^ - 1); 6 = 0.941 - 0.205(|w|°-^3^ - 1), (87) 
and e{z) = 1/1^nr(-2) - 1 = (IM - 1)(1 + zf"- 



The importance of dc and Avir arises from the fact that these quantities can be 
used to study the abundance of non hnear bound structures in the universe. 
The basic idea behind this calculation [206] is as follows: Let us consider a 
density field 5r(x) smoothed by a window function Wr of scale radius R. 
As a first approximation, we may assume that the region with 5{R,t) > 5c 
(when smoothed on the scale R at time t) will form a gravitationally bound 
object with mass M oc 'pR^ by the time t. The precise form of the M — R 
relation depends on the window function used; for a step function M = (47r/3) 
pi?^, while for a Gaussian M = (27r)^/^p-R'^- Here 6^ is a critical value for the 
density contrast given by (85) Since S{t) — D{t) for the growing mode, the 
probability for the region to form a bound structure at t is the same as the 
probability 6 > 6c[D{ti)/D{t)] at some early epoch tj. This probability can 
be easily estimated since at sufficiently early ti, the system is described by a 
Gaussian random field. This fact can be used to calculate the number density 
of bound objects leading to the result 



The quantity a here refers to the linearly extrapolated density contrast. We 
shall now describe the constraints on dark energy arising from structure for- 
mation. 



Avir(-z) = IStt^ 1 + ae\z) , 



(86) 



where 




(88) 
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5.3 Structure formation and constraints on dark energy 



Combining the initial power spectrum, P{k) — Ak'^, n ~ 1, with the transfer 
function in (75) we find that the final spectrum has the approximate form 



P{k) oc < 



Ak~^ In^ k (A; > k^ 

Ak (k < /Ceq) 



^9) 



with 27rk-^ ^ dnizeq) ^ l3{n^Yih'^y^ Mpc = 13{Th)-^h-^ Mpc [see equation 
(68)] where F = JInr^ is the shape parameter (see equation (75); we have 
assumed Qb ~ for simplicity.) From equation (89), it is clear that P{k) 
changes from an increasing function to a decreasing function at kcq , the nu- 
merical value of which is decided by the shape parameter F. Smaller values of 
Qnr and F will lead to more power at longer wavelengths. 

One of the earliest investigations which used power spectrum to determine 
was based on the APM galaxy survey [207]. This work showed that the 
existence of large scale power requires a non zero cosmological constant. This 
result was confirmed when the COBE observations fixed the amplitude of the 
power spectrum unequivocally (see section 6). It was pointed out in [208,209] 
that the COBE normalization led to a wrong shape for the power spectrum 
if we take I^nr = 1,^^a = 0, with more power at small scales than observed. 
This problem could be solved by reducing Onr and changing the shape of the 
power spectrum. Current observations favour F 0.25. In fact, an analysis of 
a host of observational data, including those mentioned above suggested [210] 
that Qaj^O even before the SN data came up. 

Another useful constraint on the models for structure formation can be ob- 
tained from the abundance of rich clusters of galaxies with masses M fa 
10^^ Mq. This mass scale corresponds to a length scale of about Sh"^ Mpc and 
hence the abundance of rich clusters is sensitive to the root-mean-square fluc- 
tuation in the density contrast at 8h~^ Mpc. It is conventional to denote this 
quantity < {Sp/pY >^/^, evaluated at Sh~^ Mpc, by ag. To be consistent with 
the observed abundance of rich clusters, equation (88) requires as ~ 0.5Qnr ^■ 
This is consistent with COBE normalization for 0.3, Q\ ~ 0.7. [Unfor- 

tunately, there is still some uncertainty about the ctq — I^nr relation. There is 
a claim [211] that recent analysis of SDSS data gives as ^ 0.33 ± 0.030^^^] 

The effect of dark energy component on the growth of linear perturbations 
changes the value of as- The results of section 5.1 translate into the fitting 
function [190] 

as - (0.50 - 0.ie)Q-^("'®), (90) 
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where e = (n-l) + (/i-0.65) and 7(0, 9) = O.21-O.22w+O.33fi + O.250. For 
constant w models with w = —1,-2/3 and —1/3, this gives cxg = 0.96,0.80 
and 0.46 respectively. Because of this effect, the abundance of clusters can 
be used to put stronger constraints on cosmology when the data for high 
redshift clusters improves. As mentioned before, linear perturbations grow 
more slowly in a universe with cosmological constant compared to the J^nr = 1 
universe. This means that chisters will be comparatively rare at high redshifts 
in a f^NR = 1 universe compared to models with cosmological constant. Only 
less than 10 per cent of massive clusters form at 2; > 0.5 in a Qnr = 1 
universe whereas almost all massive clusters would have formed by 2; f=:i 0.5 in 
a universe with cosmological constant [212,213,214,215,75]. (A simple way of 
understanding this effect is by noting that if the clusters are not in place by 
z ~ 0.5, say, they could not have formed by today in models with cosmological 
constant since there is very httle growth of fluctuation between these two 
epochs.) Hence the evolution of cluster population as a function of redshift 
can be used to discriminate between these models. 

An indirect way of measuring this abundance is through the lensing effect 
of a cluster of galaxy on extended background sources. Typically, the fore- 
ground clusters shears the light distribution of the background object and 
leads to giant arcs. Numerical simulations suggest [216] that a model with 
^^NR = 0.3, Qa = 0.7 will produce about 280 arcs which is nearly an order of 
magnitude larger than the number of arcs produced in a ^nr = 1, = 
model. (In fact, an open model with Q^b. = 0.3, = will produce about 
2400 arcs.) To use this effect, one needs a well defined data base of arcs and a 
controlled sample. At present it is not clear which model is preferred though 
this is one test which seems to prefer open model rather than a A-CDM model. 

Given the solution to (64) in the presence of dark energy, we can repeat the 
above analysis and obtain the abundance of different kinds of structures in the 
universe in the presence of dark energy. In particular this formalism can be 
used to study the abundance of weak gravitational lenses and virialized x-ray 
clusters which could act as gravitational lenses. The calculations again show 
[205] that the result is highly degenerate in w and Qnr- If ^nr is known, then 
the number count of weak lenses will be about a factor 2 smaller for w = —2/3 
compared to the ACDM model with a cosmological constant. However, if I^nr 
and w are allowed to vary in such a way that the matter power spectrum 
matches with both COBE results and abundance of x-ray clusters, then the 
predicted abundance of lenses is less than 25 per cent for —l<w< —0.4. 
It may be possible to constrain the dark energy better by comparing relative 
abundance of virialized lensing clusters with the abundance of x-ray under 
luminous lensing halos. For example, a survey covering about 50 square degrees 
of sky may be able to differentiate a ACDM model from w = —0.6 model at 
a 3(7 level. 
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The value of (Tg and cluster abundance can also be constrained via the Sunyaev- 
Zcldovich (S-Z) effect which is becoming a powerful probe of cosmological 
parameters [217]. The S-Z angular power spectrum scales as cr|(fiB^)^ and 
is almost independent of other cosmological parameters. Recently the power 
spectrum of CMBR determined by CBI and BIMA experiments (see section 6) 
showed an excess at small scales which could be interpreted as due to S-Z effect. 
If this interpretation is correct, then as{QBh/0.035)^''^^'^ = 1.04±0.12 at 95 per 
cent confidence level. This erg is on the higher side and only future observations 
can decide whether the interpretation is correct or not. The WMAP data, for 
example, leads to a more conventional value of (jg = 0.84 ± 0.04. 

Constraints on cosmological models can also arise from the modeling of damped 
Lyman-a systems [75,109,110,218,219,220] when the observational situation 
improves. At present these observations are consistent with Qnr = 0.3, a = 
0.7 model but do not exclude other models at a high significance level. 

Finally, we comment on a direct relation between S{a) and H[a). Expressing 
equation (65) in terms of H{a) will lead to the form 

a-'HH" + (3//^ + aHH')a5' = + 5) + (91) 

This can be used to determine H'^{a) from S{a) since this equation is linear 
and first order in Q{a) = H^{a) (though it is second order in S). Rewriting it 
in the form 

A(a)Q' + B{a)Q = C{a) (92) 

where 

We can integrate it to give the solution 

H\a) = iH'lflJ^^ J <;«(5^ (94) 

This shows that, given the non linear growth of perturbations S{a) as a func- 
tion of redshift and the approximate validity of spherical model, one can deter- 
mine H{a) and thus w{a) even during the nonlinear phases of the evolution. 
[A similar analysis with the linear equation (66) was done in [221], leading 
to the result which can be obtained by expanding (94) to linear order in 5.] 
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Unfortunately, this is an impractical method from observational point of view 
at present. 



6 CMBR anisotropies 

In the standard Friedmann model of the universe, neutral atomic systems form 
at a redshift of about z ^ 10^ and the photons decouple from the matter at this 
redshift. These photons, propagating freely in spacetime since then, constitute 
the CMBR observed around us today. In an ideal Friedmann universe, for a 
comoving observer, this radiation will appear to be isotropic. But if physical 
process has led to inhomogeneities in the z = 10^ spatial surface, then these 
inhomogeneities will appear as angular anisotropies in the CMBR in the sky 
today. A physical process operating at a proper length scale L on the z = 10^ 
surface will lead to an effect at an angle 9 — L/dA{z). Numerically, 

To relate the theoretical predictions to observations, it is usual to expand the 
temperature anisotropies in the sky in terms of the spherical harmonics. The 
temperature anisotropy in the sky will provide A = AT/T as a function of 
two angles 9 and ip. If we expand the temperature anisotropy distribution on 
the sky in spherical harmonics: 

AT °° 
A{9,i;) = — ^Y.^imyim{0,i^). (96) 

l,m 

all the information is now contained in the angular coefficients aim- 

If n and m are two directions in the sky with an angle a between them, the 
two-point correlation function of the temperature fluctuations in the sky can 
be defined as 

C{a) = (5(n)5(m)) ^ J2EM'm')yirn{n)Y;^,{ni). (97) 

Since the sources of temperature fluctuations are related hnearly to the density 
inhomogeneities, the coefficients aim will be random fields with some power 
spectrum. In that case < ai^aii^t > will be nonzero only iil — I' and m — ml . 
Writing 

{0'lmO'*l'm') ~ CiSli'Smm' (98) 
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and using the addition theorem of spherical harmonics, we find that 

C{a) = E ^^^^^CiPi{cosa) (99) 
I 47r 



with Ci =< \aim\'^ >. In this approach, the pattern of anisotropy is contained 
in the variation of C; with I. Roughly speaking, / oc and we can think of 
the (6,1) pair as analogue of (x, k) variables in 3-D. The Q is similar to the 
power spectrum P(k). 

In the simplest scenario, the primary anisotropics of the CMBR arise from 
three different sources, (i) The first is the gravitational potential fiuctua- 
tions at the last scattering surface (LSS) which will contribute an anisotropy 
(AT/T)| oc k^P^{k) where P^{k) oc P{k)/k'^ is the power spectrum of grav- 
itational potential 0. This anisotropy arises because photons climbing out 
of deeper gravitational wells lose more energy on the average, (ii) The sec- 
ond source is the Doppler shift of the frequency of the photons when they 
arc last scattered by moving electrons on the LSS. This is proportional to 
(AT/T)|, oc k^Py where Pv{k) oc P/A;^ is the power spectrum of the velocity 
field, (iii) Finally, we also need to take into account the intrinsic fluctuations of 
the radiation fleld on the LSS. In the case of adiabatic fluctuations, these will 
be proportional to the density fluctuations of matter on the LSS and hence 
will vary as (Ar/r)2^t oc k^P{k). Of these, the velocity fleld and the density 
fleld (leading to the Doppler anisotropy and intrinsic anisotropy described in 
(ii) and (iii) above) will oscillate at scales smaller than the Hubble radius at 
the time of decoupling since pressure support will be effective at these scales. 
At large scales, if P{k) oc k, then 

/AT\2 /AT\2 ,„ /AT 



T 



/AT\^ /AT\^ 
oc constant; j oc A;^ oc Q-"^] j oc A;^ oc 9'^ (100) 



where ^ oc A oc k^^ is the angular scale over which the anisotropy is measured. 
Obviously, the fluctuations due to gravitational potential dominate at large 
scales while the sum of intrinsic and Doppler anisotropics will dominate at 
small scales. Since the latter two are oscillatory, we sill expect an oscillatory 
behaviour in the temperature anisotropics at small angular scales. 

There is, however, one more feature which we need to take into account. The 
above analysis is valid if recombination was instantaneous; but in reality the 
thickness of the recombination epoch is about A2; ~ 80 ([222]; [44], chapter 3). 
This implies that the anisotropics will be damped at scales smaller than the 
length scale corresponding to a redshift interval of A2; = 80. The typical value 
for the peaks of the oscillation are at about 0.3 to 0.5 degrees depending on 
the details of the model. At angular scales smaller than about 0.1 degree, the 
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anisotropies are heavily damped by the thickness of the LSS. 

The fact that several different processes contribute to the structure of angular 
anisotropies make CMBR a valuable tool for extracting cosmological informa- 
tion. To begin with, the anisotropy at very large scales directly probe modes 
which are bigger than the Hubble radius at the time of decoupling and al- 
lows us to directly determine the primordial spectrum. Thus, in general, if the 
angular dependence of the spectrum at very large scales is known, one can 
work backwards and determine the initial power spectrum. If the initial power 
spectrum is assumed to be P{k) — Ak, then the observations of large angle 
anisotropy allows us to fix the amplitude A of the power spectrum [208,209]. 
Based on the results of COBE satellite [223], one finds that the amount of 
initial power per logarithmic band in k space is given by 



(This corresponds to ^4 ~ (29/i~^Mpc)^. Since the actual (AT/T) is one real- 
ization of a Gaussian random process, the observed small-/ results are subject 
to unavoidable fluctuations called the 'cosmic variance'.) This result is power- 
ful enough to rule out matter dominated, Q = 1 models when combined with 
the data on the abundance of large clusters which determines the amplitude 
of the power spectrum at R ^ Sh~^ Mpc. For example the parameter values 
h = 0.5, f^o ~ ^DM = 1,^A = 0, are ruled out by this observation when 
combined with COBE observations [208,209]. 

As we move to smaller scales we are probing the behaviour of baryonic gas cou- 
pled to the photons. The pressure support of the gas leads to modulated acous- 
tic oscillations with a characteristic wavelength at the z = 10^ surface. Regions 
of high and low baryonic density contrast will lead to anisotropies in the tem- 
perature with the same characteristic wavelength. The physics of these oscilla- 
tions has been studied in several papers in detail [224,225,226,227,228,229,230,231,232]. 
The angle subtended by the wavelength of these acoustic oscillations will lead 
to a series of peaks in the temperature anisotropy which has been detected 
[233,234]. The structure of acoustic peaks at small scales provides a reliable 
procedure for estimating the cosmological parameters. 

To illustrate this point let us consider the location of the first acoustic peak. 
Since all the Fourier components of the growing density perturbation start 
with zero amplitude at high redshift, the condition for a mode with a given 
wave vector k to reach an extremum amplitude at t = t^ec is given by 




(101) 




(102) 
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Fig. 17. The variation of the anisotropy pattern in universes with 
^NR = (0.25, 0.45, 1.0, 1.15), Oa = with the first acoustic peak moving 
from right to left. The y-axis is essentially a measure of (AT/T)^ while the x-axis 
is a measure of l/^.( Figure courtesy: S.Sethi.) 

where Cg = {dP/dpY^^ ~ (l/Vs) is the speed of sound in the baryon-photon 
fluid. At high redshifts, t{z) oc r2^p^^(l + z)~^^^ and the proper wavelength of 
the first acoustic peak scales as Apeak ~ ^dec h'^fl^K^ . The angle subtended 
by this scale in the sky depends on (Ia- If ^^nr + ^^a = 1 then the angular 
diameter distance varies as 

r^NR^ while if A = 0, it varies as ^2^^. It follows 
that the angular size of the acoustic peak varies with the matter density as 

^NR (if^A = 0), 

(103) 

Qnr' (ifl^A + ^^NR-1). 



Therefore, the angle subtended by acoustic peak is quite sensitive to ^nr if 
A = but not if f2NR+^A = 1- More detailed computations show that the mul- 

— 1/2 

tipole index corresponding to the acoustic peak scales as Ip ~ 220^^^^ if A = 
and Ip ~ 220 if JInr+^a = 1 and 0.1 < ^nr < 1. This is illustrated in figure 17 
which shows the variation in the structure of acoustic peaks when is changed 
keeping — 0. The four curves are for Q = Qnr = 0.25, 0.45, 1.0, 1.15 with 
the first acoustic peak moving from right to left. The data points on the fig- 
ures are from the first results of MAXIMA and BOOMERANG experiments 
and are included to give a feel for the error bars in current observations. It is 
obvious that the overall geometry of the universe can be easily fixed by the 
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study of CMBR anisotropy. 

The heights of acoustic peaks also contain important information. In partic- 
ular, the height of the first acoustic peak relative to the second one depends 
sensitively on Qb- However, not all cosmological parameters can be measured 
independently using CMBR data alone. For example, different models with 
the same values for {Qbm + ^a) and flBh^ will give anisotropics which are 
fairly indistinguishable. The structure of the peaks will be almost identical 
in these models. This shows that while CMBR anisotropics can, for exam- 
ple, determine the total energy density (^dm + ^a)j we will need some other 
independent cosmological observations to determine individual components. 

At present there exists several observations of the small scale anisotropics in 
the CMBR from the balloon flights, BOOMERANG [233], MAXIMA [234], 
and from radio telescopes CBI [235], VSA [236], DASI [237,238] and — most 
recently — from WMAP [239] . These CMBR data has been extensively ana- 
lyzed [239,59,60,223,235,236,240,241,242,243,244,245] in isolation as well as in 
combination with other results. (The information about structure formation 
arises mainly from galaxy surveys like SSRS2, CfA2 [246], ICRS [247], Abell- 
ACO cluster survey [248], IRAS-PSC z [243], 2-D survey [249,242] and the 
Sloan survey [250] .) While there is some amount of variations in the results, 
by and large, they support the following conclusions. 

• The data strongly supports a /c = model of the universe [245] with Qtot = 
1.00±[]:g^ from the pre-MAP data and Qtot = 1-02 ± 0.02 from the WMAP 
data. 

• The CMBR data before WMAP, when combined with large scale structure 
data, suggest Qnr = 0.29 ± 0.05 ± 0.04 [59,60,245,251]. The WMAP result 
[239] is consistent with this giving 0.27 ±0.04. The initial power spectrum is 
consistent with being scale invariant and the pre-MAP value is n = 1.02 ± 
0.06 ± 0.05 [59,60,245]. The WMAP gives the spectral index at A; = 0.05 
Mpc~^ to be 0.93 ±0.03. In fact, combining 2dF survey results with CMBR 
suggest [252] 0.7 independent of the supernova results. 

• A similar analysis based on BOOMERANG data leads to Qtot — 1-02 ± 
0.06(scc for example, [241]). Combining this result with the HST constraint 
[49] on the Hubble constant h = 0.72 ± 0.08, galaxy clustering data as well 
SN observations one gets Ha = 0.62%%, = 0.55%%^ and Ha = 0-7^-l%7 
respectively [253]. The WMAP data gives h = 0.7llg':g|. 

• The analysis also gives an independent handle on baryonic density in the 
universe which is consistent with the BBN value: The pre-MAP result was 
^Bh'^ = 0.022 ± 0.003 [59,60]. (This is gratifying since the initial data had 
an error and gave too high a value [254].) The WMAP data gives flsh^ — 
0.0224 ± 0.0009. 
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There has been some amount of work on the effect of dark energy on the CMBR 
anisotropy [255,256,257,258,259,260,261,262,263]. The shape of the CMB spec- 
trum is relatively insensitive to the dark energy and the main effect is to alter 
the angular diameter distance to the last scattering surface and thus the posi- 
tion of the first acoustic peak. Several studies have attempted to put a bound 
on w using the CMB observations. Depending on the assumptions which were 
invoked, they aU lead to a bound broadly in the range of w < —0.6. (The pre- 
liminary analysis of WMAP data in combination with other astronomical data 
sets suggest w < —0.78 at 95 per cent confidence limit.) At present it is not 
clear whether CMBR anisotropics can be of significant help in discriminating 
between different dark energy models. 



7 Reinterpreting the cosmological constant 

It is possible to attack the cosmological constant problem from various other 
directions in which the mathematical structure of equation (3) is reinterpreted 
differently. Though none of these ideas have been developed into a successful 
formal theory, they might contain ingredients which may eventually provide a 
solution to this problem. Based on this hope, we shall provide a brief descrip- 
tion of some of these ideas. (In addition to these ideas, there is extensive lit- 
erature on several different paradigms for attacking the cosmological constant 
problem based on: (i) Quantum field theory in curved spacetime [264,265,266], 
(n) quantum cosmological considerations [267], (in) models of infiation [268], 
(iv) string theory inspired ideas [269], and (v) effect of phase transitions [270].) 

7.1 Cosmological constant as a Lagrange multiplier 

The action principle for gravity in the presence of a cosmological constant 

can be thought of as a variational principle extremizing the integral over R, 
subject to the condition that the 4-volume of the universe remains constant. 
To implement the constraint that the 4-volume is a constant, one will add a 
Lagrange multiplier term which is identical in structure to the second term in 
the above equation. Hence, mathematically, one can think of the cosmological 
constant as a Lagrange multiplier ensuring the constancy of the 4-volume of 
the universe when the metric is varied. 
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If we take this interpretation seriously, then it is necessary to specify the 4- 
vohime of the universe before the variation is performed and determine the 
cosmological constant so that the 4- volume has this specified volume. A Fried- 
mann model with positive cosmological constant in Minkowski space will lead 
to a finite 3-volume proportional to A~^/^ on spatial integration. (To achieve 
this, we should use the coordinates in which the spatial sections are closed 
3-spheres.) The time integration, however, has an arbitrary range and one 
needs to restrict the integration to part of this range by invoking some phys- 
ical principle. If we take this to be typically the age of the universe, then 
we will obtain a time dependent cosmological constant A{t) with A(t)H(t)~'^ 
remaining of order unity. While this appears to be a conceptually attractive 
idea, it is not easy to implement it in a theoretical model. In particular, it is 
difficult to obtain this as a part of a generally covariant theory incorporating 
gravity. 

7.2 Cosmological constant as a constant of integration 

Several people have suggested modifying the basic structure of general relativ- 
ity so that the cosmological constant will appear as a constant of integration. 
This does not solve the problem in the sense that it still leaves its value un- 
determined. But this changes the perspective and allows one to think of the 
cosmological constant as a non dynamical entity [271,272]. 

One simple way of achieving this is to assume that the determinant g of gab 
is not dynamical and admit only those variations which obeys the condition 
g'^^^Qah = in the action principle. This is equivalent to eliminating the trace 
part of Einstein's equations. Instead of the standard result, we will now be led 
to the equation 

i?^ - \51R = SttG (tI - i^r) (105) 

which is just the traceless part of Einstein's equation. The general covari- 
ance of the action, however, implies that T.f' = and the Bianchi identi- 
ties [Rl — ^6lR).i = continue to hold. These two conditions imply that 
diR = —SiiGdiT requiring R + SuGT to be a constant. Calling this constant 
(—4 A) and combining with equation (105), we get 

Ri - ^5iR - 5iA = SnCTi (106) 

which is precisely Einstein's equation in the presence of cosmological constant. 
In this approach, the cosmological constant has nothing to do with any term 
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in the action or vacuTim fluctuations and is merely an integration constant. 
Like any other integration constant its value can be fixed by invoking suitable 
boundary conditions for the solutions. 

There are two key difficulties in this approach. The first, of course, is that it 
still does not give us any handle on the value of the cosmological constant and 
all the difficulties mentioned earlier still exists. This problem would have been 
somewhat less serious if the cosmological constant was strictly zero; the pres- 
ence of a small positive cosmological constant makes the choice of integration 
constant fairly arbitrary. The second problem is in interpreting the condition 
that g must remain constant when the variation is performed. It is not easy to 
incorporate this into the logical structure of the theory. (For some attempts 
in this direction, see [273].) 



7.3 Cosmological constant as a stochastic variable 

Current cosmological observations can be interpreted as showing that the ef- 
fective value value of A (which will pick up contributions from all vacuum 
energy densities of matter fields) has been reduced from the natural value of 
L~p^ to L~^[LpHqY where Hq is the current value of the Hubble constant. One 
possible way of thinking about this issue is the following [274]: Let us assume 
that the quantum micro structure of spacetime at Planck scale is capable of 
readjusting itself, soaking up any vacuum energy density which is introduced 
— like a sponge soaking up water. If this process is fully deterministic and 
exact, all vacuum energy densities will cease to have macroscopic gravitational 
effects. However, since this process is inherently quantum gravitational, it is 
subject to quantum fluctuations at Planck scales. Hence, a tiny part of the 
vacuum energy will survive the process and will lead to observable effects. One 
may conjecture that the cosmological constant we measure corresponds to this 
small residual fluctuation which will depend on the volume of the spacetime 
region that is probed. It is small, in the sense that it has been reduced from 
Lp^ to Lp^^LpHoY, which indicates the fact that fluctuations — when mea- 
sured over a large volume — is small compared to the bulk value. It is the 
wetness of the sponge we notice, not the water content inside. 

This is particularly relevant in the context of standard discussions of the con- 
tribution of zero-point energies to cosmological constant. The correct theory is 
likely to regularise the divergences and make the zero point energy finite and 
about Lp^. This contribution is most likely to modify the microscopic struc- 
ture of spacetime (e.g if the spacetime is naively thought of as due to stacking 
of Planck scale volumes, this will modify the stacking or shapes of the volume 
elements) and will not affect the bulk gravitational field when measured at 
scales coarse grained over sizes much bigger than the Planck scales. 
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Given a large 4-volume V of the spacetime, we will divide it into M cTibes of 

size (Ax)^ and label the cubes by n = 1,2, , M. The contribution to the 

path integral amplitude A, describing long wavelength limit of conventional 
Einstein gravity, can be expressed in the form 



^ = [exp(ci(i?L^ 



) + ■•■) 



exp —J- I a xy 



-g{RLl) (107) 



where we have indicated the standard continuum limit. (In conventional units 
Ci = (IGtt)^^.) Let us now ask how one could modify this result to describe 
the ability of spacetime micro structure to readjust itself and absorb vacuum 
energy densities. This would require some additional dynamical degree of free- 
dom that will appear in the path integral amplitude and survive in the classical 
limit. It can be shown that [274] the simplest implementation of this feature 
is by modifying the standard path integral amplitude [exp(ci(i?Lp) + ■■■)] 
by a factor [0(a;„)/0o] where (j){x) is a scalar degree of freedom and 0o is a 
pure number introduced to keep this factor dimensionless. In other words, we 
modify the path integral amplitude to the form: 



modify 



n 



0(a^n 
00 



(108) 



In the long wavelength limit, the extra factor in (108) will lead to a term of 
the form 



IJexp 



■ exp 



ijp 



I 



d^Xy 




(109) 



Thus, the net effect of our assumption is to introduce a 'scalar field potential' 
y (0) = -Lp^ In (0/0o) in the semi classical limit. It is obvious that the rescal- 
ing of such a scalar field by — > g0 is equivalent to adding a cosmological 
constant with vacuum energy —Lj,'^ In q. Alternatively, any vacuum energy can 
be re absorbed by such a rescaling. The fact that the scalar degree of freedom 
occurs as a potential in (109) without a corresponding kinetic energy term 
shows that its dynamics is unconventional and non classical. 



The above description in terms of macroscopic scalar degree of freedom can, of 
course, be only approximate. Treated as a vestige of a quantum gravitational 
degrees of freedom, the cancellation cannot be precise because of fluctua- 
tions in the elementary spacetime volumes. These fluctuations will reappear 
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as a "small" cosmological constant because of two key ingredients: (i) discrete 
spacetime structure at Planck length and (ii) quantum gravitational uncer- 
tainty principle. 

To show this, we use the fact noted earlier in section 7.1 that the net cosmo- 
logical constant can be thought of as a Lagrange multiplier for proper volume 
of spacetime in the action functional for gravity. In any quantum cosmologi- 
cal models which leads to large volumes for the universe, phase of the wave 
function will pick up a factor of the form 



from (104), where V is the four volume. Treating (Aggf/STrLp, V) as conjugate 
variables {q,p), we can invoke the standard uncertainty principle to predict 
A A Ki 87rLp/AV. Now we use the ear her assumption regarding the microscopic 
structure of the spacetime: Assume that there is a zero point length of the 
order of Lp so that the volume of the imiverse is made of a large number 
{N) of cells, each of volume (aLp)'^ where a is a numerical constant. Then 
V = N{aLpY, implying a Poisson fluctuation AV ~ \/V(aI/p)^ and leading 
to 



This will give Qa = (Stt/Sq;^) which will — for example — lead to — (2/3) 
if q; = 2y^. Thus Planck length cutoff (UV limit) and volume of the universe 
(IR limit) combine to give the correct AA. (A similar result was obtained 
earlier in [273] based on a different model.) The key idea, in this approach, is 
that A is a stochastic variable with a zero mean and ffuctuations. It is the rms 
fluctuation which is being observed in the cosmological context. 

This has three implications: First, FRW equations now need to be solved 
with a stochastic term on the right hand side and one should check whether 
the observations can still be explained. The second feature is that stochastic 
properties of A need to be described by a quantum cosmological model. If the 
quantum state of the universe is expanded in terms of the eigenstates of some 
suitable operator (which does not commute the total four volume operator), 
then one should be able to characterize the ffuctuations in each of these states. 
Third, and most important, the idea of a cosmological constant arising as a 
fluctuation makes sense only if the bulk value is rescaled away. 

The non triviality of this result becomes clear when we compare it with few 
other alternative ways of estimating the ffuctuations — none of which gives 
the correct result. The first alternative approach is based on the assumption 




(110) 




(111) 
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that one can associate an entropy S = {Ah/'^L']:>) with compact space time 
horizons of area (Wc will discuss this idea in detail in section 10). A popular 
interpretation of this result is that horizon areas are quantized in units of L|> 
so that S is proportional to the number of bits of information contained in 
the horizon area. In this approach, horizon areas can be expressed in the form 
Ah = ApN where Ap (x Lp is a quantum of area and N is an integer. Then the 
fluctuations in the area will be AAh = Ap\^ = y/ApAn- Taking Ah oc A^^ 
for the de Sitter horizon, we find that A A oc H^{HLp) which is a lot smaller 
than what one needs. Further, taking Ah oc rjj, we find that Arn oc Lp; that 
is, this result essentially arises from the idea that the radius of the horizon is 
uncertain within one Planck length. This is quite true, of course, but does not 
lead to large enough fluctuations. 

A more sophisticated way of getting this (wrong) result is to relate the fluc- 
tuations in the cosmological constant to that of the volume of the universe is 

by using a canonical ensemble description for universes of proper Euclidean 
4- volume [275]. Writing V = V/87rLp and treating V and A as the relevant 
variables, one can write a partition function for the 4-volume as 

oo 

Z{V) = J g{A)e-^^dA (112) 





Taking the analogy with standard statistical mechanics (with the correspon- 
dence V ^ P and A —>£■), we can evaluate the fluctuations in the cosmological 
constant in exactly the same way as energy fluctuations in canonical ensemble. 
(This is done in several standard text books; see, for example, [276] p. 194.) 
This will give 

(AA?-^- C- --V^^^ (113) 

^^^^ - T^' ^ - W/v) ~ ^ dV ^^^^^ 



where C is the analogue of the spcciflc heat. Taking the 4-vohimc of the 
universe to be V = = 9bA^^ where 6 is a numerical factor and using 

V = (V/SnLl) we get A oc LpV-V2. it follows from (113) that 

(AA)^ = ^ = ^-^{LpHr (114) 



In other words A A oc H^{HLp), which is the same result from area quantiza- 
tion and is a lot smaller than the cosmologically signiflcant value. 

Interestingly enough, one could do slightly better by assuming that the horizon 
radius is quantized in units of Planck length, so that th = = NLp. 

This will lead to the fluctuations Arn — V'>^hLp and using rn — H'^ oc 
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A~^/^, we get AA oc H'^{HLpY/'^ — larger than (114) but still inadequate. 
In summary, the existence of two length scales H^^ and Lp allows different 
results for AA depending on how exactly the fluctuations are characterized 
(Ay oc -s/TV, oc s/N or Ar^ oc \/N). Hence the result obtained above in 
(111) is non trivial. 

These conclusions stress, among other things, the difference between fluctua- 
tions and the mean values. For, if one assumes that every patch of the universe 
with size Lp contained an energy £'p, then a universe with characteristic size 
H"^ will contain the energy E = {Ep/Lp)H~^. The corresponding energy 
density will be pa = {E/H~^) = {H/Lp)^ which leads to the correct result. 
But, of course, we do not know why every length scale Lp should contain an 
energy Ep and — more importantly — contribute coherently to give the total 
energy. 

7.4 Anthropic interpretation of the cosmological constant 

The anthropic principle [277,278] is an interpretational paradigm which ar- 
gues that, while discussing the origin of physical phenomena and the values of 
constants of nature, we must recognize the fact that only certain combination 
and range of values will lead to the existence of intelligence observers in the 
universe who could ask questions related to these issues. This paradigm has 
no predictive power in the sense that none of the values of the cosmological 
parameters were ever predicted by this method. ^ In fact some cosmologists 
have advocated the model with I^nr = 1, ^^a = strongly and later — when 
observations indicated Qa 7^ 1 — have advocated the anthropic interpretation 
of cosmological constant with equal fluency. This is defended by the argument 
that not all guiding principles in science (Darwinian evolution, Plate tectonics, 
....) need to be predictive in order to be useful. In this view point, anthropic 
principle is a back drop for discussing admittedly complicated conceptual is- 
sues. Within this paradigm there have been many attempts to explain (after 
the fact) the values of several fundamental constants with varying degree of 
success. 

In the context of cosmological constant, the anthropic interpretation works 
as follows. It is assumed that widely disparate values for the constants of na- 
ture can occur in an ensemble of universes (or possibly in different regions 
of the universe causally unconnected with each other). Some of these values 

^ Some advocates of the anthropic principle cite Fred Hoyle predicting the existence 
of excited state of carbon nucleus, thereby leading to efficient triple alpha reaction 
in stellar nucleosynthesis, as an example of a prediction from anthropic principle; it 
is very doubtful whether Hoyle applied anthropic considerations in arriving at this 
conclusion. 
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for constants of nature — and in particular for the cosmological constant — 
will lead broadly to the kind of universe we seem to live in. This is usually 
characterized by formation of: (i) structures by gravitational instability, (ii) 
stars which act as gravitationally bound nuclear reactors that synthesize the 
elements and distribute them and (iii) reasonably complex molecular struc- 
tures which could form the basis for some kind of life form. Showing that such 
a scenario can exist only for a particular range of values for the cosmological 
constant is considered an explanation for the value of cosmological constant 
by the advocates of anthropic principle. (More sophisticated versions of this 
principle exist; see, for example [279], and references cited therein.) 

The simplest constraint on the cosmological constant is that it should not be 
so high as to cause rapid expansion of the universe early on preventing the 
formation of galaxies [280]. If the energy density of the cosmological constant 
has to be less than that of energy density of matter at the redshift 2;gai(~ 4) 
at which galaxy formation takes place, then we must have 

< (1 + ^gal)' ~ 125 (115) 



This gives a bound on Qa which is "only" a couple of orders of magnitude 
larger than what is observed. 

More formally, one could ask: What is the most probable value of if it is 
interpreted as the value that would have been observed by the largest number 
of observers [281,282]? Since a universe with Qa ~ ^nr will have more galaxies 
than one with a universe with ^a ~ IO^^^nR) one could argue that most 
observers will measure a value I^a ~ ^nr- The actual probability dP for 
measuring a particular value for fl\ in the range (^a, QA + d^A) is the product 
{dP/dflA) — Q{QA)V{i^A) where V is the a priori probability measure for a 
specific value of f^A in a member of an ensemble of universes (or in a region 
of the universe) and Q{Qa) is the average number of galaxies which form in a 
universe with a given value of Oa- There has been several attempts to estimate 
these quantities (see, for example, [283,284]) but all of them are necessarily 
speculative. The first — and the most serious — difficulty with this approach is 
the fact that we simply do not have any reliable way of estimating V; in fact, if 
we really had a way of calculating it from a fundamental theory, such a theory 
probably would have provided a deeper insight into the cosmological constant 
problem itself. The second issue has to do with the dependence of the results 
on other parameters which describe the cosmological structure formation (like 
for example, the spectrum of initial perturbations). To estimate Q one needs 
to work in a multi parameter space and marginalize over other parameters 
— which would involve more assumptions regarding the priors. And finally, 
anthropic paradigm itself is suspect in any scientific discussion, for reasons 
mentioned earlier. 
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7.5 Probabilistic interpretation of the cosmological constant 

It is also possible to produce more complex scenarios which could justify the 
small or zero value of cosmological constant. One such idea, which enjoyed 
popularity for a few years [285,286,287,288], is based on the conjecture that 
quantum wormholes can change the effective value of the observed constants 
of nature. The wave function of the universe, obtained by a path integral over 
all possible spacetime metrics with wormholes, will receive dominant contri- 
butions from those configurations for which the effective values of the physical 
constants extremize the action. Under some assumptions related to Euclidean 
quantum gravity, one could argue that the configurations with zero cosmolog- 
ical constant will occur at late times. It is, however, unlikely that the assump- 
tions of Euclidean quantum gravity has any real validity and hence this idea 
must be considered as lacking in concrete justification. 



8 RelcLxation mechanisms for the cosmological constant 

One possible way of obtaining a small, non-zero, cosmological constant at the 
present epoch of the universe is to make the cosmological constant evolve in 
time due to some physical process. At a phcnomcnological level this can be 
done either by just postulating such a variation and explore its consequences or 
— in a slightly more respectable way — by postulating a scalar field potential 
as described in section 4. These models, however, cannot explain why a bare 
cosmological constant [the first term on the right hand side of (7)] is zero. To 
tackle this issue, one can invoke some field [usually a scalar field] which directly 
couples to the cosmological constant and decreases its "effective value". We 
shall now examine two such models. 

The key idea is to introduce a field which couples to the trace T = of 
the energy momentum tensor. If T depends on and vanishes at some value 
(p = (po, then (p will evolve towards (p = (po ai which T = 0. This equilibrium 
solution will have zero cosmological constant[289,290,291,292]. While this idea 
sounds attractive, there are general arguments as to why it does not work in 
the simplest context [4] . 

A related attempt was made by several authors, [289,293,294,295], who cou- 
pled the scalar field directly to it! which, of course, is proportional to T because 
of Einstein's equations. Generically, these models have the Lagrangian 



L^\^^d,cPd^cP + 



IGttG 



1 



{R - 2A) - U{(p)R 



(116) 
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The field equations of this model has flat spacctimc solutions at (f) = (po pro- 
vided U{(f)o) = oo. Unfortunately, the effective gravitational constant in this 
model evolves as 

G 

= 1 + 167rGC/(0o) ^^^^^ 



and vanishes as U — > oo. Hence these models are not viable. 

The difficulty in these models arise because they do not exphcitly couple the 
trace of the Tab of the scalar field itself. Handling this consistently [296] leads 
to a somewhat different model which we will briefly describe because of its 
conceptual interest. 

Consider a system consisting of the gravitational fields gab, radiation fields, and 
a scalar field which couples to the trace of the energy-momentum tensor of 
all fields, including its own. The zeroth order action for this system is given 

by 

^(0) = A,,,, + Af + ASI + Aradn (118) 



where 



Agrav = {16ttG)-^ J R^^d'^x - J A^^d^x, (119) 
<t>%y/^9 d'x; 42 = V J Tf{cj>/<j>^),/^ d'x (120) 



Here, we have explicitly included the cosmological constant term and 77 is a 
dimcusioiilcss number which 'switches on' the interaction. In the zeroth order 
action, T represents the trace of all fields other than (f). Since the radiation 
field is traceless, the only zeroth-order contribution to T comes from the A 
term, so that we have T = 4A. The coupling to the trace is through a function 
/ of the scalar field, and one can consider various possibilities for this function. 
The constant 0o converts to a dimensionless variable, and is introduced for 
dimensional convenience. 

To take into account the back-reaction of the scalar field on itself, we must 
add to T the contribution — —(f)^(f)i of the scalar field. If we now add to T 
in the interaction term A^l further modifies T^^ This again changes T^. Thus 

to arrive at the correct action an infinite iteration will have to be performed 
and the complete action can be obtained by summing up all the terms. (For a 
demonstration of this iteration procedure, see [297,298].) The full action can 
be found more simply by a consistency argument. 
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Since the effect of the iteration is to modify the expression for and ^a, we 
consider the following ansatz for the full action: 



(121) 



Here q;(0) and are functions of to be determined by the consistency 
requirement that they represent the effect of the iteration of the interaction 
term. (Since radiation makes no contribution to T, we expect A^^d to remain 
unchanged.) The energy-momentum tensor for and A is now given by 



(122) 



so that the total trace is Ttot = 4q;(0)A — I3{(t))(t)^(j)i. The functions q;(^) and 
/3(0) can now be determined by the consistency requirement 



- j a{ct>)K^d^x+]^ J (3{(l))(l>'(l)i^d^x 

= - 1 AV^ d'^ + lj 'P''P^^9 d^x + r]J T,,J{(j)/(j)o)V^d^x{123) 
Using Ttot and comparing terms in the above equation we find that 

a(0) = [1 + Arif]-', (5{4>) = [1 + 2r;/]-^ (124) 

Thus the complete action can be written as 



1 

+ 7^ 



#4: 



2 J l + 2n/ 



+ 4nf 
^d'^x + A^^ 



(125) 



(The same action would have been obtained if one uses the iteration proce- 
dure.) The action in (125) leads to the following field equations, 



/3(0) u> 



A 



SttG 



(126) 
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(127) 



Here, □ stands for a covariant d'Lambertian, J'tra<:eiess jg ^j^^ stress tensor of 
all fields with traceless stress tensor and a prime denotes differentiation witfi 
respect to (f). 

In the cosmological context, this reduces to 



'^+T'^^^'^rT^ + "2^ (1 + W 



+ A: 



57r 



1 0' A 1 Po' 

21 + 2?7/ 87rG(l + 4?7/) 



(128) 
(129) 



It is obvious that the effective cosmological constant can decrease if / increases 
in an expanding universe. The result can be easily generalized for a scalar field 
with a potential by replacing A by V{(j)). This model is conceptually attractive 
since it correctly accounts for the coupling of the scalar field with the trace of 
the stress tensor. 

The trouble with this model is two fold: (a) If one uses natiual initial conditions 
and do not fine tune the parameters, then one does not get a viable model, 
(b) Since the scalar field couples to the trace of all sources, it also couples to 
dust-like matter and "kills" it, making the universe radiation dominated at 
present. This reduces the age of the universe and could also create difficulties 
for structure formation. These problems can be circumvented by invoking a 
suitable potential K(0) within this model [299]. However, such an approach 
takes away the naturalness of the model to certain extent. 



9 Geometrical structure of the de Sitter spacetime 



The most symmetric vacuum solution to Einstein's equation, of course, is the 
flat spacetime. If we now add the cosmological constant as the only source of 
curvature in Einstein's equation, the resulting spacetime is also highly sym- 
metric and has an interesting geometrical structure. In the case of a positive 
cosmological constant, this is the de Sitter manifold and in the case of negative 
cosmological constant, it is known as anti-de Sitter manifold. We shall now 
discuss some features of the former, corresponding to the positive cosmological 
constant. (For a nice, detailed, review of the classical geometry of de Sitter 
spacetime, see [300].) 
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To understand the geometrical structure of the de Sitter spacetime, let us 
begin by noting that a spacetime with the source T^" = Pa^6 must have 3- 
dimensional section which are homogeneous and isotropic. This will lead us 
to the Einstein's equations for a FRW universe with cosmological constant as 
source 



a' k SttG 



This equation can be solved with any of the following three forms of {k, a{t)) 
pair. The first pair is the spatially flat universe with (k — 0,a — e^*). The 
second corresponds to spatially open universe with (A; = — 1, a = sinh Ht) 
and the third will be (A; = +1, a = cosh Ht). Of these, the last pair gives 
a coordinate system which covers the full de Sitter manifold. In fact, this 
is the metric on a 4-dimensional hyperboloid, embedded in a 5 dimensional 
Minkowski space with the metric 

ds"^ = dt^ - dx^ - dy^ - dz^ - dv^, (131) 
The equation of the hyperboloid in 5-D space is 

i2_^2_^2_^2_^2^_^-2_ (^32) 



We can introduce a parametric representation of the hyperbola with the four 
variables (r, x, ^, 0) where 



x = H ^ cosh(i7r) sin X sin 6* cos 0; y = H ^ cosh (i^r) sin x sin 6* sin 0; 
z = H^^ cosh(iJr) sin x cos 9; 



V — H ^ cosh(i/r) cosx; 
t = H-^smh{HT). 



(133) 



This set, of course, satisfies (132). Using (131), we can compute the metric 
induced on the hyperboloid which — when expressed in terms of the four 
coordinates {T,x,9,(f)) — is given by 

ds^ = dr^ - H-^ cosh^(Hr) [dx^ + sin^ x(d9^ + sin^ ^d0^)] , (134) 



This is precisely the de Sitter manifold with closed spatial sections. 

All the three forms of FRW universes with A; = 0, ±1 arise by taking different 
cuts in this 4-dimensional hyperboloid embedded in the 5-dimensional space- 
time. Since two of these dimensions (corresponding to the polar angles 9 and 
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(j)) merely go for a ride, it is more convenient (for visualization) to work with 
a 3-dimensional spacetime having the metric 

ds'^ = - dx^ - dv^. (135) 

instead of the 5-dimcnsional metric (131). Every point in this 3-dimcnsional 
space corresponds to a 2-sphere whose coordinates 6 and are suppressed for 
simplicity. The (1 + 1) de Sitter spacetime is the 2-dimensional hyperboloid 
[instead of the four dimensional hyperboloid of (132)] with the equation 

-x"^ -v'^^ -H-'^ (136) 

embedded in the 3-dimensional space with metric (135). The three different 
coordinate systems which are natural on this hyperboloid are the following: 

• Closed Spatial Sections: This is obtained by introducing the coordinates 
t = smh.{HT);x = cosh(i7r) sinx; v = cosh(i7r) cosx on the 
hyperboloid, in terms of which the induced metric on the hyperboloid has 
the form 

ds^ = dr^ - H-^ cosh\HT)dx^ (137) 

This is the two dimensional de Sitter space which is analogous to the 4- 
dimensional case described by (134). 

• Open Spatial Sections: These are obtained by using the coordinates t = 

sinh(i?r) cosh^; x = sinh(i?T) sinh^; v = cosh(i7r) on the 
hyperboloid in terms of which the induced metric on the hyperboloid has 
the form 

ds'^ = dr"^ - H-^ smh^ {HT)di^ (138) 

• Flat Spatial Sections: This corresponds to the choice t — sinh(i7T) + 
leading to the metric 

ds'^ = dr^ - exp(2i7r)d{2 (^^39) 

This covers one half of the de Sitter hyperboloid bounded by the null rays 
t + x = Q. 

All these metrics have an apparent time dependence. But, in the absence of 

any source other than cosmological constant, there is no preferred notion of 
time and the spacetime manifold cannot have any intrinsic time dependence. 
This is indeed true, in spite of the expansion factor a{t) ostensibly depending 
on time. The translation along the time direction merely slides the point on the 
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surface of the hyperboloid. [This is obvious in the coordinates (/c = 0, a oc e^*) 
in which the time translation t ^ t + e merely rescales the coordinates by 
(exp He).] 

The time independence of the metric can be made explicit in another set of 
coordinates called 'static coordinates'. To motivate these coordinates, let us 
note that a spacetime with only cosmological constant as the source is certainly 
static and possesses spherical symmetry. Hence we can also express the metric 
in the form 

ds^ = e'^dt^ - e^dr^ - r\de^ + sin^ edcj)^) (140) 

where p and A are functions of r. The Einstein's equations for this metric has 
the solution e'^ — = (1 — H'^r'^) leading to 

ds^ - (1 - H^r^)dt^ - , ^''l^ - r\dQ'' + sin^ Qd^^) (141) 

(1 — H^r^) 

This form of the metric makes the static nature apparent. This metric also 
describes a hyperboloid embedded in a higher dimensional flat space. For 
example, in the (1 + 1) case (with 9^ suppressed) this metric can be obtained 
by the following parameterization of the hyperboloid in equation (136): 

t = {H-'^ - r^)^/^ sinh(i/T); v = (//"^ - r'^f'^ cosh(HT); x^r (142) 

The key feature of the manifold, revealed by equation (141) is the existence 
of a horizon at r = H~^. It also shows that i is a time-like coordinate only in 
the region r < H~^. 

The structure of the metric is very similar to the Schwarzschild metric: 

ds'' = (l - ^) dt' - - r\de^ + sin^ ^#2) (143) 

Both the metrics (143) and (141) are spherically symmetric with g^Q = ^{l/gu) 
Just as the Schwarzschild metric has a horizon at r = 2M (indicated by 
5*00 0,5111 — > 00), the de Sitter metric also has a horizon at r = H~^. 
Prom the slope of the hght cones (dt/dr) = ±(1 — if^r^)~^ [corresponding 
to ds — — d6 — d(f) in (142)] it is clear that signals sent from the region 
r < cannot go beyond the surface r = H~^. 

This feature, of course, is independent of the coordinate system used. To see 
how the horizon in de Sitter universe arises in the FRW coordinates, let us 
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recall the equation governing the propagation of light signals between the 
events (ii,ri) and {t,r): 



[ ^ (144) 

J VI - kx-" J a{t') ^ ' 



Consider a photon emitted by an observer at the origin at the present epoch 
(ri = 0, ^1 = to). The maximum coordinate distance xh reached by this photon 
as t ^ oo is determined by the equation 



" dx r dt' 



i «(^')' 



(145) 



If the integral on the right hand side diverges as t ^ oo, then, in the same 
limit, Xff ^ oo and an observer can send signals to any event provided (s)he 
waits for a sufficiently long time. But if the integral on the right hand side 
converges to a finite value as t — > oo, then there is a finite horizon radius 
beyond which the observer's signals will not reach even if (s)he waits for infinite 
time. In the de Sitter universe with A; = and a{t) = e^*, xh = H~^e~^^'^; the 
corresponding maximum proper distance up to which the signals can reach is 
rn — a{to)xH — H~^. Thus we get the same result in any other coordinate 
system. 

Since the result depends essentially on the behaviour of a{t) as t — > oo, it will 
persist even in the case of a universe containing both non relativistic matter 
and cosmological constant. For example, in our universe, we can ask what is 
the highest redshift source from which we can ever receive a light signal, if 
the signal was sent today. To compute this explicitly, consider a model with 
^NR + f^A = 1. Let us assume that light from an event at {th-iZh) reaches 
r = at z = giving 

to ZH 

f dt f dz ,^ 

^H= / ^ = / TTTT- 146 

I a{t) I Ho[l- Qnr + J1nr(1 + zff' 



If we take th to be the size of the horizon, then it also follows that the light 
emitted today from this event will just reach us at i = oo. This gives 

oo , , 

dt r dz 



[^=[ ^ (147) 

J a{t) I Ho [1 - l^NR + ^^nr(1 + ^)31'/' 



to 



Equating the two expressions, we get an implicit expression for zh- If ^^nr ~ 
0.3, the limiting redshift is quite small: zh ~ 1.8. This imphes that sources 
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with z > zh can never be influenced by light signals from us in a model with 
cosmological constant [301,302]. 



10 Horizons, temperature and entropy 

In the description of standard cosmology appears as a parameter like, say, 
the Hubble constant Hq. There is, however, a signiflcant difference between 
these two parameters as far as fundamental physics is concerned. The exact 
numerical value of h is not of major concern to fundamental physics. But, the 
non-zero value for Q,\ signifies the existence of an exotic form of energy density 
with negative pressure which is a result of deep significance to the whole of 
physics. We shall now take up an important aspect of the cosmological constant 
which is somewhat different in spirit compared to the results covered so far. 

It turns out that the universe with a non-zero value for cosmological constant 
behaves in many ways in a manner similar to a black hole. Just as the black 
hole has close links with thermodynamics (like having a finite temperature, 
entropy etc.) the dc Sitter universe also possesses thermodynamic features 
which makes it peculiar and important in understanding the cosmological 
constant. This thermodynamic relationship of the cosmological constanthas 
not been adequetely explored or integrated into the standard cosmological 
description so far. But since it is likely to have a important implications for 
the eventual resolution of the cosmological constantproblem, we shall provide 
a fairly self contained description of the same. 

One of the remarkable features of classical gravity is that it can wrap up regions 
of spacetime thereby producing surfaces which act as one way membranes. 
The classic example is that of Schwarzschild black hole of mass M which has 
a compact spherical surface of radius r = 2M that act as a horizon. Since 
the horizon can hide information — and information is deeply connected with 
entropy — one would expect a fundamental relationship between gravity and 
thermodynamics. [There is extensive literature in this subject and our citation 
will be representative rather than exhaustive; for a text book discussion and 
earlier references, sec [303]; for a recent review, see [304].] As we saw in the 
last section, the de Sitter universe also has a horizon which suggests that de 
Sitter spacetime will have non trivial thermodynamic features [305] . 

This result can be demonstrated mathematically in many different ways of 
which the simplest procedure is based on the relationship between temperature 
and the Euclidean extension of the spacetime. To see this connection, let 
us recall that the mean value of some dynamical variable f{q) in quantum 
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statistical mechanics can be expressed in the form 

<f>-^12j rE{Q)fiQ)MQ)e-^'' dq (148) 

where (psiQ) is the stationary state eigen function of the Hamihonian with 
= E(f)E,P = i^/T) is the inverse temperature and Z{P) is the partition 
function. This expression calculates the mean value < E\f\E > in a given 
energy state and then averages over a Boltzmann distribution of energy states 
with the weightage exp{—f3E). On the other hand, the quantum mechani- 
cal kernel giving the probability amplitude for the system to go from the state 
q at time t — to the state q' at time t is given by 

K{q', t; q,0)^Y: 4>l{q')M<l)e-'''' (149) 

E 

Comparing (148) and (149) we find that the thermal average in (148) can be 
obtained by 

<f>=^jdq Kiq, -il3- q, 0)f{q) (150) 

in which we have done the following: (i) The time coordinate has been analyt- 
ically continued to imaginary values with it = r. (ii) The system is assumed 
to exhibit periodicity in the imaginary time r with period (3 in the sense that 
the state variable q has the same values at r = and at r = /3. These con- 
siderations continue to hold even for a field theory with q denoting the field 
configuration at a given time. If the system, in particular the Greens functions 
describing the dynamics, are periodic with a period p in imaginary time, then 
one can attribute a temperature T — (1/p) to the system. It may be noted 
that the partition function Z{P) can also be expressed in the form 

^(/3) = Ee"^'' = JdQK{q, -iP;q,0)^ | P?exp[-A^(g, /3; g, 0)] (151) 

The first equality is the standard definition for Z{j3)\ the second equality fol- 
lows from (149) and the normalization of (Pe{q)] the last equality arises from 
the standard path integral expression for the kernel in the Euclidean sector 
(with Ae being the Euclidean action) and imposing the periodic boundary 
conditions. (It is assumed that the path integral measure Vq includes an in- 
tegration over q.) We shall have occasion to use this result later. Equations 
(150) and (151) represent the relation between the periodicity in Euclidean 
time and temperature. 
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Spacetimcs with horizons possess a natural analytic continuation from Minkowski 
signature to the Euclidean signature with t ^ t = it. \i the metric is peri- 
odic in r, then one can associate a natural notion of a temperature to such 
spacetimes. For example, the de Sitter manifold with the metric (134) can be 
continued to imaginary time arriving at the metric 



- ds^ = dr^ + H-^ cos^ Ht 



dx^ + sin^ xid^'^ + sin^ 9d4>'^ 



(152) 



which is clearly periodic in r with the period {27i/H). [The original metric was 
a 4-hyperboloid in the 5-dimensional space while equation (152) represents a 
4-sphere in the 5-dimensional space.] It follows that de Sitter spacetime has a 
natural notion of temperature T — {H/2ti) associated with it. 

It is instructive to see how this periodicity arises in the static form of the 
metric in (141). Consider a metric of the form 

ds'' = f{r)de -^-dLl (153) 

where dL^ denotes the transverse 2-dimensional metric and /(r) has a simple 
zero at r = th- Near r = th, we can expand /(r) in a Taylor series and obtain 
/(r) ~ B{r — th) where B = f'{rH)- The structure of the metric in (153) 
shows that there is a horizon aXr = th- Further, since the general relativistic 
metric reduces to Qqq ~ (1 -|- 207v) in the Newtonian limit, where is the 
Newtonian gravitational potential, the quantity 

^ = WN{rH)\ - \\9UrH)\ - \\nrH)\ = \\B\ (154) 

can be interpreted as the gravitational attraction on the surface of the horizon 
— usually called the surface gravity. Using the form /(r) 2/t(r — th) near 
the horizon and shifting to the coordinate ^ = [2«;~^(r — r^)]^^^ the metric 
near the horizon becomes 

ds" ^ K^i^de - di^ - dLi (155) 



The Euclidean continuation t ^ t — it now leads to the metric 

- rfs^ ^ ^^d{KTf + d^'^ + dLl (156) 



which is essentially the metric in the polar coordinates in the t — ^ plane. For 
this metric to be well defined near the origin, kt should behave like an angular 
coordinate 9 with periodicity 27r. Therefore, we require all well defined physical 
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quantities defined in this spacetime to have a periodicity in r with the period 
(27r/|K|). Thus, all metrics of the form in (153) with a simple zero for /(r) 
leads to a horizon with temperature T = |fi;|/27r = |/'(rij)|/47r. In the case 
of de Sitter spacetime, this gives T — (i?/27r); for the Schwarzschild metric, 
the corresponding analysis gives the well known temperature T — (I/SttM) 
where M is the mass of the black-hole. 



10.1 The connection between thermodynamics and spacetime geometry 



The existence of one-way membranes, however, is not necessarily a feature 
of gravity or curved spacetime and can be induced even in flat Minkowski 
spacetime. It is possible to introduce coordinate charts in Minkowski space- 
time such that regions are separated by horizons, a familiar example being 
the coordinate system used by a uniformly accelerated frame (Rindler frame) 
which has a non-compact horizon. The natural coordinate system (t, x, y, z) 
used by an observer moving with a uniform acceleration g along the x-axis is 
related to the inertial coordinates (T, X, Y, Z) by 

gT = yr+2^sinh(^t); (1 + gX) = ^ I + 2gx cosh{gt)- (157) 



and Y — y; Z — z. The metric in the accelerated frame will be 

ds^^{l + 2gx)dt^--^j^^^-dy^-dz^ (158) 



which has the same form as the metric in (153) with f[x) = (1 -|- 2gx). This 
has a horizon at x — —l/2g with the surface gravity k = g and temperature 
T — {g/27r). All the horizons are imphcitly defined with respect to certain class 
of observers; for example, a suicidal observer plunging into the Schwarzschild 
black hole will describe the physics very differently from an observer at in- 
finity. From this point of view, which we shall adopt, there is no need to 
distinguish between observer dependent and observer independent horizons. 
This allows a powerful way of describing the thermodynamical behaviour of 
all these spacetimes (Schwarzschild, de Sitter, Rindler ....) at one go. 

The Schwarzschild, de Sitter and Rindler metrics are symmetric under time re- 
versal and there exists a 'natural' definition of a time symmetric vacuum state 
in all these cases. Such a vacuum state will appear to be described a thermal 
density matrix in a subregion TZ of spacetime with the horizon as a boundary. 
The QFT based on such a state will be manifestedly time symmetric and will 
describe an isolated system in thermal equilibrium in the subregion TZ. No 
time asymmetric phenomena like evaporation, outgoing radiation, irreversible 
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changes etc can take place in this situation. We shall now describe how this 
arises. 

Consider a (D+1) dimensional flat Lorentzian manifold S with the signature 
(+, — , — , ...) and Cartesian coordinates where A — (0, 1, 2, D). A four 
dimensional sub-manifold V in this (D+1) dimensional space can be defined 
through a mapping = Z^i^x"") where with a = (0, 1,2,3) are the four 
dimensional coordinates on the surface. The flat Lorentzian metric in the 
(D+1) dimensional space induces a metric gabi^"") on the four dimensional 
space which — for a wide variety of the mappings Z^ — Z^{x"') — will 
have the signature (+, — , — , — ) and will represent, in general, a curved four 
geometry. The quantum theory of a free scalar fleld in S is well deflned in 
terms of the, say, plane wave modes which satisfy the wave equation in S. A 
subset of these modes, which does not depend on the 'transverse' directions, 
will satisfy the corresponding wave equation in V and will depend only on x". 
These modes induce a natural QFT in V. We are interested in the mappings 
Z^ = Z'^i^x"') which leads to a horizon in V so that we can investigate the 
QFT in spacetimes with horizons using the free, flat spacetime, QFT in S 
([309] [304]). 

For this purpose, let us restrict attention to a class of surfaces defined by the 
mappings Z"^ = Z^i^x"") which ensures the following properties for V: (i) The 
induced metric Qat has the signature (+, — , — , — ). (ii) The induced metric is 
static in the sense that goa = and all QabS are independent of x^. [The Greek 
indices run over 1,2,3.] (iii) Under the transformation — > x°±i(7r/g'), where 
g is a non zero, positive constant, the mapping of the coordinates changes as 
ZO ^ -z°, Z^ -Z^ and Z^ Z^ for A = 2,...,D.lt will turn out that the 
four dimensional manifolds defined by such mappings possess a horizon and 
most of the interesting features of the thermodynamics related to the horizon 
can be obtained from the above characterization. Let us first determine the 
nature of the mapping = Z^(x") = Z'^{t, x) such that the above conditions 
are satisfied. 

The condition (iii) above singles out the spatial coordinate Z^ from the others. 
To satisfy this condition we can take the mapping Z"^ — Z'^{t, r,6,(j)) to be 
of the form Z° = Z°{t,r),Z^ = Z\t,r),Z^ = Z^{r,9,(f)) where Z^ denotes 
the transverse coordinates Z"^ with A = {2,...,D). To impose the condition 
(ii) above, one can make use of the fact that S possesses invariance under 
translations, rotations and Lorentz boosts, which are characterized by the 
existence of a set of = (1/2)(I^ + 1){D + 2) KiUing vector fields C^(Z^). 
Consider any linear combination of these Killing vector fields which is time 
like in a region of S. The integral curves to this vector field V"^ will define time 
like curves in S. If one treats these curves as the trajectories of a hypothetical 
observer, then one can set up the proper Fermi- Walker transported coordinate 
system for this observer. Since the four velocity of the observer is along the 
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Killing vector field, it is obvious that the metric in this coordinate system will 
be static [310]. In particular, there exists a Killing vector which corresponds 
to Lorentz boosts along the direction that can be interpreted as rotation in 
imaginary time coordinate allowing a natural realization of (iii) above. Using 
the property of Lorentz boosts, it is easy to see that the transformations of the 
form Z'^ = lf{rY^'^ sinhty't; — ±//(r)^/^ cosh (^t will satisfy both conditions 
(ii) and (iii) where (/,(?) are constants introduced for dimensional reasons 
and /(r) is a given function. This map covers only the two quadrants with 
\Z^\ > \Z^\ with positive sign for the right quadrant and negative sign for the 
left. To cover the entire Z^) plane, we will use the full set 

Z^ = lf{rY''^smhgt] Z^ ^ ±lf{rfl'^ cosh gt (for \Z^\>\Z^\) (159) 
Z°-±/[-/(r)]^/2cosh^t; Z^ = l[-f{r)Y'^sixi\igt (for \Z^\ < |Z°|) 

The inverse transformations corresponding to (159) are 

ff{r) = {Z'f-{Zy; gt^t^nh-'iZ'/Z') (160) 

Clearly, to cover the entire two dimensional plane of — oo < {Z^, Z^) < +oo, 
it is necessary to have both /(r) > and /(r) < 0. The pair of points (Z°, Z^) 
and {—Z^, —Z^) are mapped to the same {t, r) making this a 2-to-l mapping. 
The null surface Z^ — ±Z^ is mapped to the surface /(r) — 0. 

The transformations given above with any arbitrary mapping for the trans- 
verse coordinate Z-^ — Z^{r, 9, 4>) will give rise to an induced metric on T> of 
the form 

ds'' = f{r){lgfde - J dr^ - dLl (161) 

where dL\ depends on the form of the mapping Z-^ — Z-^{r, 9, 0). This form 
of the metric is valid in all the quadrants even though we will continue to work 
in the right quadrant and will comment on the behaviour in other quadrants 
only when necessary. It is obvious that the in general, is curved and has a 
horizon at /(r) = 0. 

As a specific example, let us consider the case of (D+l)=6 with the coordinates 
(ZO, Z\ Z^, Z^, Z^ Z^) = (ZO, Z^ Z^, i?, e, $) and consider a mapping to 4- 
dimensional subspace in which: (i) The (Z°, Z^) are mapped to (i, r) as before; 
(ii) the spherical coordinates (i?, ©, $) in S are mapped to standard spherical 
polar coordinates in V: (r, 9, Lp) and (iii) we take Z^ to be an arbitrary function 
of r: Z^ = Qit^)- This leads to the metric 

ds' = A{r)dt' - B{r)dr' - r'dQ',,; (162) 
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with 

A{r)^{lgrf; S(r) = 1 + g'^ + (163) 



Equation (162) is the form of a general, sphericaUy symmetric, static metric in 
4-dimension with two arbitrary functions /(r), q{r). Given any specific metric 
with A{r) and B{r), equations (163) can be solved to determine /(r), q{r). As 
an example, let us consider the Schwarzschild solution for which we will take 
/ = 4 [1 — (l/r)]; the condition goo — (l/gu) now determines q{r) through the 
equation 




(164) 




Though the integral cannot be expressed in terms of elementary functions, it 
is obvious that g(r) is well behaved everywhere inchiding at r = /. The trans- 
formations {Z^,Z^) {t,r);Z'^ q{i^)', {Z'"' , Z'^ , Z^) {r,6,Lp) thus provide 
the embedding of Schwarzschild metric in a 6-dimensional space. [This result 
was originally obtained by Frondsal [311]; but the derivation in that paper is 
somewhat obscure and does not bring out the generality of the situation] . As 
a corollary, we may note that this procedure leads to a spherically symmetric 
Schwarzschild-like metric in arbitrary dimension, with the 2-sphere in (162) 
replaced any A'"-sphere. 

The choice Ig — l,/(r) = [1 — (r//)^] will provide an embedding of the de 
Sitter spacetime in 6-dimensional space with Z^ — r, (Z^, Z^, Z^) (r, 6, 0). 
Of course, in this case, one of the coordinates is actually redundant and — 
as we have seen earlier — one can achieve the embedding in a 5-dimensional 
space. A still more trivial case is that of Rindler metric which can be obtained 
with D=3, Ig — 1, /(r) — 1 + 2gr; in this case, the "embedding" is just a 
reparametrization within four dimensional spacetime and — in this case — r 
runs in the range (—00,00). The key point is that the metric in (161) is fairly 
generic and can describe a host of spacetimes with horizons located at / = 0. 
We shall discuss several features related to the thermodynamics of the horizon 
in the next few sections. 



84 



10.2 Temperature of horizons 

There exists a natural definition of QFT in the original (D + 1) -dimensional 
space; in particular, we can define a vacuum state for the quantum field on 
the = surface, which coincides with the t = surface. By restricting the 
field modes (or the field configurations in the Schrodingcr picture) to depend 
only on the coordinates in T>, we will obtain a quantum field theory in T> 
in the sense that these modes will satisfy the relevant field equation defined 
in v. In general, this is a complicated problem and it is not easy to have a 
choice of modes in S which will lead to a natural set of modes in V. We can, 
however, take advantage of the arguments given in the last section — that all 
the interesting physics arises from the (Z°, Z^) plane and the other transverse 
dimensions are irrelevant near the horizon. In particular, solutions to the wave 
equation in <S which depends only on the coordinates Z^ and Z^ will satisfy 
the wave equation in V and will depend only on {t, r). Such modes will define 
a natural s— wave QFT in V. The positive frequency modes of the above 
kind (varying as exp{—iQZ^) with a; > 0.) will be a specific superposition 
of negative (varying as e*'^*) and positive (varying as e~*'^*) frequency modes 
in T> leading to a temperature T — {g/27r) in the 4-dimensional subspace 
on one side of the horizon. There are several ways of proving this result, all 
of which depend essentially on the property that under the transformation 
t — >• t ± (in/g) the two coordinates Z^ and Z^ reverses sign. 

Consider a positive frequency mode of the form Fq{Z^,Z^) oc exp[—iQZ'^ + 
iPZ^] with r2 > 0. These set of modes can be used to expand the quantum 
field thereby defining the creation and annihilation operators An, A^: 

ct>{Z\ Z') = YX^^MZ\ Z^) + Af*^{Z\ Z')] (166) 
n 

The vacuum state defined by A^jvac >= corresponds to a globally time sym- 
metric state which will be interpreted as a no particle state by observers using 
Z^ as the time coordinate. Let us now consider the same mode which can be 
described in terms of the (i, r) coordinates. Being a scalar, this mode can be ex- 
pressed in the 4-dimcnsional sector in the form Fn(t, r) = FQ[Z^{t, r), Z^(t, r)\. 
The Fourier transform of FQ,{t, r) with respect to t will be:, 

oo 

Ka{u,r)^ I dte+'^'Fn[Z''{t,r),Z\t,r)]; (-oo < u; < oo) (167) 

— oo 

Thus a positive frequency mode in the higher dimension can only be expressed 
as an integral over uj with uj ranging over both positive and negative values. 
However, using the fact that t ^ t — {iTi/g) leads to Z^ — > — Z^ — > —Z^, 
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it is easy to show that 



(168) 



This allows us to write the inverse relation to (167) as 



-oo 
oo 



I ^ [i^a(a;,r)e— * + e--/^ira(a;,r)e^-* 



(169) 



The term with represents the contribution of negative frequency modes in 
the the 4-D spacetime to the pure positive frequency mode in the embedding 
spacetime. A field mode of the embedding spacetime containing creation and 
annihilation operators {Aq, Aq) can now be represented in terms of the cre- 
ation and annihilation operators (a(^,aj,) appropriate to the {t,r) coordinates 
as 



AnFn + A^F^ = [{a^ + ^e"""/^) Kne~'^' + h.c.' 



duo 1 



a^Kae-''^' + h.c. 



(170) 



where A^^; is a normalization constant. Identifying = N^{Aq, + e^'^'^/^A^) 
and using the conditions [at^,a^] = 1, [Aq,^]^] = 1 etc., we get = [1 — 
exp(— 27ra;/g')]~-'^/^. It follows that the number of a— particles in the vacuum 
defined by An|vac >= is given by 



< vacla^OijIvac 



(171) 



This is a Planckian spectrum with temperature T = 5'/27r. The key role in the 
derivation is played by equation (168) which, in turn, arises from the analytic 
properties of the spacetime under Euclidean continuation. 



10.3 Entropy and energy of de Sitter spacetime 



The best studied spacetimes with horizons are the black hole spacetimes. (For 
a sample of references, see [312,313,314,315,316,317,318,319,320,321]). In the 
simplest context of a Schwarzschild black hole of mass M, one can attribute an 
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energy E = M, temperature T — (SttM)"^ and entropy S = {1/4:){Ah/Lp) 
where Ah is the area of the horizon and Lp = [Gh/c'Y^'^ is the Planck length. 
(Hereafter, we will use units with G = h = c = 1.) These are clearly related 
by the thermodynamic identity TdS — dE, usually called the first law of 
black hole dynamics. This result has been obtained in much more general 
contexts and has been investigated from many different points of view in the 
literature. The simplicity of the result depends on the following features: (a) 
The Schwarzschild metric is a vacuum solution with no pressure so that there 
is no PdV term in the first law of thermodynamics, (b) The metric has only 
one parameter M so that changes in all physical parameters can be related to 
dM. (c) Most importantly, there exists a well defined notion of energy E to 
the spacetime and the changes in the energy dE can be interpreted in terms of 
the physical process of the black hole evaporation. The idea can be generalized 
to other black hole spacetimes in a rather simple manner only because of well 
defined notions of energy, angular momentum etc. 

Can one generalize the thermodynamics of horizons to cases other than black 
holes in a straight forward way ? In spite of years of research in this field, 
this generalization remains non trivial and challenging when the conditions 
listed above are not satisfied. To see the importance of the above conditions, 
we only need to contrast the situation in Schwarzschild spacetime with that 
of de Sitter spacetime: 

• As we saw in section 10.2, the notion of temperature is well defined in the 
case of de Sitter spacetime and we have T — H/2t^ where is the radius 
of the de Sitter horizon. But the correspondence probably ends there. A 
study of literature shows that there exist very few concrete calculations 
of energy, entropy and laws of horizon dynamics in the case of de Sitter 
spacetimes, in sharp contrast to BH space times. 

• There have been several attempts in the hterature to define the concept of 
energy using local or quasi-local concepts (for a sample of references, see 
[322,323,324,325,326,327,328,329,330]). The problem is that not aU defini- 
tions of energy agree with each other and not all of them can be applied to 
de Sitter type universes. 

• Even when a notion of energy can be defined, it is not clear how to write and 
interpret an equation analogous to dS = [dE/T) in this spacetime, espe- 
cially since the physical basis for dE would require a notion of evaporation 
of the de Sitter universe. 

• Further, we know that de Sitter spacetime is a solution to Einstein's equa- 
tions with a source having non zero pressure. Hence one would very much 
doubt whether TdS is indeed equal to dE. It would be necessary to add a 
PdV term for consistency. 

All these suggest that to make any progress, one might require a local approach 
by which one can define the notion of entropy and energy for spacetimes with 
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horizons. This conclusion is strengthened further by the following argument: 
Consider a class of spherically symmetric spacetimes of the form 

ds^ = f{r)dt' - fir)-^dr^ - r^dO'' + sin' ed(t>^) (172) 

If /(r) has a simple zero at r = a with /'(a) = B remaining finite, then this 
spacetime has a horizon at r = a. Spacetimes like Schwarzschild or de Sitter 
have only one free parameter in the metric (like M or H~^) and hence the 
scaling of all other thermodynamical parameters is uniquely fixed by purely 
dimensional considerations. But, for a general metric of the form in (172), with 
an arbitrary /(r), the area of the horizon (and hence the entropy) is deter- 
mined by the location of the zero of the function /(r) while the temperature 
— obtained from the periodicity considerations — is determined by the value 
of /'(r) at the zero. For a general function, of course, there will be no relation 
between the location of the zero and the slope of the function at that point. It 
will, therefore, be incredible if there exists any a priori relationship between 
the temperature (determined by /' ) and the entropy (determined by the zero 
of /) even in the context of horizons in spherically symmetric spacetimes. If 
wc take the entropy to be S = tto,' (where /(a) = determines the radius of 
the horizon) and the temperature to be T = |/'(a)|/47r (determined by the 
periodicity of Euclidean time), the quantity TdS — {l/2)\f'{a)\ada will de- 
pend both on the slope f'{a) as well as the radius of the horizon. This implies 
that any local interpretation of thermodynamics will be quite non trivial. 

Finally, the need for local description of thermodynamics of horizons becomes 
crucial in the case of spacetimes with multiple horizons. The strongest and the 
most robust result we have, regarding spacetimes with a horizon, is the notion 
of temperature associated with them. This, in turn, depends on the study of 
the periodicity of the Euclidean time coordinate. This approach does not work 
very well if the spacetime has more than one horizon like, for example, in the 
Schwarzschild-de Sitter metric which has the form in (172) with 

/(r) = (l - ^ - ^V) (173) 
This spacetime has two horizons at r± with 

^+"V3 — 3 — ' ^ 

where cos a; = —?)^/?)MH~^ . (The parameter x is in the range (vr, (3/2)7r] 
and we assume that < 27A'PH~'^ < 1.) Close to either horizon the space- 
time can be approximated as Rindler. Since the surface gravities on the two 
horizons are different, we get two different temperatures T± — |/'(r±)|/47r. 




(174) 
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To maintain invariance under it ^ it + (3 (with some finite (3) it is nec- 
essary that j3 is an integer multiple of both 47r/|/'(r+)| and 47r/|/'(r_)| so 
that P = (47rn-|-/|/'(r±)|) where n± are integers. Hence the ratio of surface 
gravities |/'(r+)|/|/'(r_)| = {rij^/nJ) must be a rational number. Though irra- 
tionals can be approximated by rationals, such a condition definitely excludes 
a class of values for M if i7 is specified and vice versa. It is not clear why the 
existence of a cosmological constant should imply something for the masses of 
black holes (or vice versa). Since there is no physical basis for such a condi- 
tion, it seems reasonable to conclude that these difficulties arise because of our 
demanding the existence of a finite periodicity (5 in the Euclidean time coordi- 
nate. This demand is related to an expectation of thermal equilibrium which 
is violated in spacetimes with multiple horizons having different temperatures. 

If even the simple notion of temperature falls apart in the presence of multiple 
horizons, it is not likely that the notion of energy or entropy can be defined 
by global considerations. On the other hand, it will be equally strange if we 
cannot attribute a temperature to a black hole formed in some region of the 
universe just because the universe at the largest scales is described by a de 
Sitter spacetime, say. One is again led to searching for a local description of 
the thermodynamics of all types of horizons. We shall now see how this can 
be done. 

Given the notion of temperature, there are two very different ways of defin- 
ing the entropy: (1) In statistical mechanics, the partition function Z[I3) of 
the canonical ensemble of systems with constant temperature is related 
to the entropy S and energy E by Z^P) oc exp(S' — PE). (2) In classical 
thermodynamics, on the other hand, it is the change in the entropy, which 
can be operationally defined via dS = dE/T{E). Integrating this equation 
will lead to the function S{E) except for an additive constant which needs to 
be determined from additional considerations. Proving the equality of these 
two concepts was nontrivial and — historically — led to the unification of 
thermodynamics with mechanics. 

In the case of time symmetric state, there will be no change of entropy dS 
and the thermodynamic route is blocked. It is, however, possible to construct a 
canonical ensemble of a class of spacetimes and evaluate the partition function 
Z{P). For spherically symmetric spacetimes with a horizon at r = /, the parti- 
tion function has the generic form Z oc cxp[S — f3E], where S = (l/4)47r/^ and 
\E\ = (//2). This analysis reproduces the conventional result for the black hole 
spacetimes and provides a simple and consistent interpretation of entropy and 
energy for de Sitter spacetime, with the latter being given by = —(1/2)H~^. 
In fact, it is possible to write Einstein's equations for a spherically symmetric 
spacetime as a thermodynamic identity TdS — dE = PdV with T, S and E 
determined as above and the PdV term arising from the source [308]. We shall 
now discuss some of these issues. 
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Consider a class of spacetimes with the metric 

ds^ = f{r)df - f{r)-Hr^ - dLl 



(175) 



where /(r) vanishes at some surface r = I, say, with /'(/) = B remaining 
finite. When dL]^ — r'^dS^ with [0 < r < oo], equation (175) covers a variety 
of spherically symmetric spacetimes with a compact horizon at r = Z. Since 
the metric is static, Euclidean continuation is trivially effected by t — > r = it 
and an examination of the conical singularity near r = a [where /(r) k, B{r — 
a)] shows that r should be interpreted as periodic with period {3 = 47r/|i?| 
corresponding to the temperature T — \B\/At^. Let us consider a set S of such 
metrics in (175) with the restriction that [/(a) = 0,f'{a) — B] but /(r) is 
otherwise arbitrary and has no zeros. The partition function for this set of 
metrics S is given by the path integral sum 

Z{P) = Y.^M-M9)) = Eexp \-T^ fdrj d^x^RE[f{T)]\ (176) 

geS geS I iOTT 7 7 



where Einstein action has been continued in the Euclidean sector and we have 
imposed the periodicity in r with period P = 4:7i/\B\. The sum is restricted 
to the set S of all metrics of the form in (175) with the behaviour [/(a) — 
0, f'{a) — B] and the Euchdean Lagrangian is a functional of /(r). The spatial 
integration will be restricted to a region bounded by the 2-spheres r — a 
and r = b, where the choice of b is arbitrary except for the requirement 
that within the region of integration the Lorentzian metric must have the 
proper signature with t being a time coordinate. The remarkable feature is 
the form of the Euclidean action for this class of spacetimes. Using the result 
R = Vlf — {2/r'^){d/dr) [r(l — /)] vahd for metrics of the form in (175), a 
straight forward calculation shows that 



" / * [-[rVr + 2[r(l - /)]'] = '^{a'B - 2a] + /'(6)](177) 



where Q depends on the behaviour of the metric near r = b and we have 
used the conditions [/(a) = 0, /'(a) = B]. The sum in (176) now reduces 
to summing over the values of [fib),f'{b)] with a suitable (but unknown) 
measure. This sum, however, will only lead to a factor which we can ignore in 
deciding about the dependence of Z{f3) on the form of the metric near r = a. 
Using P = An/B (and taking B > 0, for the moment) the final result can be 
written in a very suggestive form: 



Z(/3) = Zoexp 



oc exp [S{a) - PE{a)] 



(178) 
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with the identifications for the entropy and energy being given by: 



S = i(4.a=) = \a 



( 



hori: 



) 



1/2 



E= -a = 
2 



■izon 



(179) 



horizon j 



IGtt 



In the case of the Schwarzschild black hole with a = 2M, the energy turns 
out to he E — (a/2) = M which is as expected. (More generally, E — 
(^horizon/167r)^/^ corrcspouds to the so called 'irreducible mass' in BH space- 
times [331].) Of course, the identifications S = (AttM^), E = M,T^ (I/SttM) 
are consistent with the result dE = TdS in this particular case. 

The above analysis also provides an interpretation of entropy and energy in the 
case of de Sitter universe. In this case, /(r) = (1 — i7^r^), a = H~^, B = —2H. 
Since the region where t is time like is "inside" the horizon, the integral for 
ylfi in (177) should be taken from some arbitrary value r = b to r = a with 
a > 6. So the horizon contributes in the upper limit of the integral introducing 
a change of sign in (177). Further, since B < 0, there is another negative sign 
in the area term from PB oc B/\B\. Taking all these into account we get, in 
this case. 



giving S = (l/4)(47ra2) = (l/4)Aorizon and E = -(1/2) H'K These def- 



initions do satisfy the relation TdS — PdV — dE when it is noted that 
the dc Sitter universe has a non zero pressure P = —p\ = —E/V asso- 
ciated with the cosmological constant. In fact, if we use the "reasonable" 
assumptions S = (l/4)(47riJ"2), oc H'^ and E = -PV in the equa- 
tion TdS — PdV — dE and treat E as an unknown function of H, we get 
the equation H'^{dE/dH) = —{2>EH + 1) which integrates to give precisely 
E = —{1/2)H^^. (Note that we only needed the proportionality, V oc in 
this argument since PdV oc (dV/V). The ambiguity between the coordinate 
and proper volume is irrelevant.) 

A peculiar feature of the metrics in (175) is worth stressing. This metric will 
satisfy Einstein's equations provided the the source stress tensor has the form 



Tj: = T^ = (e(r)/87r); = Tf = (/x(r)/87r). The Einstein's equations now 



The remarkable feature about the metric in (175) is that the Einstein's equa- 
tions become linear in /(r) so that solutions for different e(r) can be super- 




reduce to: 



(181) 
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posed. Given any e(r) the solution becomes 

V 

a 1 



/(r) = l---- / eirydr (182) 
r r J 



with a being an integration constant and ^{r) is fixed by e(r) through: yu(r) = 
e+(l/2)re'(r). The integration constant a in (182) is chosen such that /(r) = 
at r = a so that this surface is a horizon. Let us now assume that the solution 
(182) is such that f{r) — Q aX r — a with /'(a) = B finite leading to leading 
to a notion of temperature with (3 — (47r/|S|). Prom the first of the equations 
(181) evaluated at r = a, we get 

\Ba-]-^-\e{a)a' (183) 



It is possible to provide an interesting interpretation of this equation which 
throws light on the notion of entropy and energy. Multiplying the above equa- 
tion by da and using e = SttT^, it is trivial to rewrite equation (183) in the 
form 




= -T;(a)[47raVa (184) 



Let us first consider the case in which a particular horizon has f (a) = B > 
so that the temperature is T = B/An. Since f{a) = 0, /'(a) > 0, it follows that 
/ > for r > a and / < for r < a; that is, the "normal region" in which t is 
time like is outside the horizon as in the case of, for example, the Schwarzschild 
metric. The first term in the left hand side of (184) clearly has the form of 
TdS since we have an independent identification of temperature from the 
periodicity argument in the local Rindler coordinates. Since the pressure is 
P = —T^, the right hand side has the structure of PdV or — more relevantly 
— is the product of the radial pressure times the transverse area times the 
radial displacement. This is important because, for the metrics in the form 
(175), the proper transverse area is just that of a 2-sphere though the proper 
volumes and coordinate volumes differ. In the case of horizons with B = 
f'{a) > which we arc considering (with da > 0), the volume of the region 
where / < will increase and the volume of the region where / > will 
decrease. Since the entropy is due to the existence of an inaccessible region, 
dV must refer to the change in the volume of the inaccessible region where 
/ < 0. We can now identify T in TdS and P in PdV without any difficulty and 
interpret the remaining term (second term in the left hand side) as dE = da/2. 
We thus get the expressions for the entropy S and energy E (when i? > 0) to 
be the same as in (179). 
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Using (184), we can again provide an interpretation of entropy and energy in 
the case of de Sitter universe. In this case, /(r) = (1 — H'^r'^), a = H~^, B — 
—2H < so that the temperature — which should be positive — is T = 
|/'(a)|/(47r) = {-B)/An. For horizons with B = /'(a) < (hke the de Sitter 
horizon) which we are now considering, /(a) = 0, /'(a) < 0, and it follows 
that / > for r < a and / < for r > a; that is, the "normal region" in 
which t is time like is inside the horizon as in the case of, for example, the de 
Sitter metric. Multiplying equation (184) by (—1), we get 

-B ,fl . o\ 1 . . . /47r 



47r 



d (^47ra=^) + = Tl{a)d (yo^) = P{-dV) (185) 



The first term on the left hand side is again of the form TdS (with positive 
temperature and entropy) . The term on the right hand side has the correct sign 
since the inaccessible region (where / < 0) is now outside the horizon and the 
volume of this region changes by {—dV). Once again, we can use (185) to iden- 
tify [308] the entropy and the energy: S = (l/4)(47ra^) = ( 1/4) ^horizon ; = 
—{1/2)H^^. These results agree with the previous analysis. 



10.4 Conceptual issues in de Sitter thermodynamics 



The analysis in the last few sections was based on a strictly static 4-dimensional 
spacetime. The black hole metric, for example, corresponds to an eternal black 
hole and the vacuum state which we constructed in section 10.2 corresponds 
to the Hartle-Hawking vacuum [332] of the Schwarzschild spacetime, describ- 
ing a black hole in thermal equilibrium. There is no net radiation flowing to 
inflnity and the entropy and temperature obtained in the previous sections 
were based on equilibrium considerations. 

As we said before, there are two different ways of defining the entropy. In 
statistical mechanics, the entropy S{E) is related to the degrees of freedom [or 
phase volume] g{E) by S{E) = lng[E). Maximization of the phase volume for 
systems which can exchange energy will then lead to equality of the quantity 
T{E) = (dS/dE)"^ for the systems. It is usual to identify this variable as 
the thermodynamic temperature. The analysis of BH temperature based on 
Hartle-Hawking state is analogous to this approach. 

In classical thermodynamics, on the other hand, it is the change in the entropy 
which can be operationally defined via dS — dE/T{E). Integrating this equa- 
tion will lead to the function S{E) except for an additive constant which needs 
to be determined from additional considerations. This suggests an alternative 
point of view regarding thermodynamics of horizons. The Schwarzschild met- 
ric, for example, can be thought of as an asymptotic limit of a metric arising 
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from the collapse of a body forming a black-hole. While developing the QFT 
in such a spacetime containing a collapsing black-hole, we need not maintain 
time reversal invariance for the vacuum state and — in fact — it is more 
natural to choose a state with purely in-going modes at early times like the 
Unruh vacuum state [333]. The study of QFT in such a spacetime shows that, 
at late times, there will exist an outgoing thermal radiation of particles which 
is totally independent of the details of the collapse. The temperature in this 
case will be T(M) = I/SttM, which is the same as the one found in the 
case of the state of thermal equilibrium around an "eternal" black-hole. In 
the Schwarzschild spacetime, which is asymptotically flat, it is also possible 
to associate an energy E = M with the black-hole. Though the calculation 
was done in a metric with a fixed value of energy E = M, it seems reason- 
able to assume that — as the energy flows to inflnity at late times — the 
mass of the black hole will decrease. If we make this assumption — that the 
evaporation of black hole will lead to a decrease of M — then one can inte- 
grate the equation dS = dM/T{M) to obtain the entropy of the black- hole 
to be S" = AttM^ = {1/A){A/Ll) where A = ATi{2Mf is the area of the event 
horizon and Lp = {Gh/c^Y^'^ is the Planck length.^ The procedure outlined 
above is similar in spirit to the approach of classical thermodynamics rather 
than statistical mechanics. 

Once it is realized that only the asymptotic form of the metric matters, we can 
simplify the above analysis by just choosing a time asymmetric vacuum and 
working with the asymptotic form of the metric with the understanding that 
the asymptotic form arose due to a time asymmetric process (like gravitational 
collapse). In the case of black hole spacetimes this is accomplished — for 
example — by choosing the Unruh vacuum [333]. The question arises as to 
how our unifled approach fares in handling such a situation which is not time 
symmetric and the horizon forms only asymptotically as i — > oo. 

There exist analogues for the collapsing black-hole in the case of de Sitter 
(and even Rindler) [308] . The analogue in the case of dc Sitter spacetime 
will be an FRW universe which behaves like a de Sitter universe only at late 
times [like in equation (27); this is indeed the metric describing our universe if 
= 0.7, Qnr = 0.3]. Mathematically, we only need to take a{t) to be a func- 
tion which has the asymptotic form exp{Ht) at late times. Such a spacetime 
is, in general, time asymmetric and one can choose a vacuum state at early 



^ This integration can determine the entropy only up to an additive constant. To 
fix this constant, one can make the additional assumption that S should vanish 
when M = 0. One may think that this assumption is eminently reasonable since the 
Schwarzschild metric reduces to the Lorcntzian metric when M ^ 0. But note that 
in the same limit of M 0, the temperature of the black-hole diverges !. Treated 
as a limit of Schwarzschild spacetime, normal flat spacetime has inflnite — rather 
than zero — temperature. 
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times in such a way that a thermal spectrum of particles exists at late times. 
Emboldened by the analogy with black-hole spacetimes, one can also directly 
construct quantum states (similar to Unruh vacuum of black-hole spacetimes) 
which are time asymmetric, even in the exact de Sitter spacetime, with the 
understanding that the de Sitter universe came about at late times through a 
time asymmetric evolution. 

The analogy also works for Rindler spacetime. The coordinate system for an 
observer with time dependent acceleration will generalize the standard Rindler 
spacetime in a time dependent manner. In particular, one can have an observer 
who was inertial (or at rest) at early times and is uniformly accelerating at late 
times. In this case an event horizon forms at late times exactly in analogy with 
a collapsing black-hole. It is now possible to choose quantum states which are 
analogous to the Unruh vacuum - which will correspond to an inertial vacuum 
state at early times and will appear as a thermal state at late times. The 
study of different 'vacuum' states shows [308] that radiative flux exists in the 
quantum states which are time asymmetric analogues of the Unruh vacuum 
state. 

A formal analysis of this problem will involve setting up the in and out vacua 

of the theory, evolving the modes from t = —oo to t = +00, and computing 
the Bogoliubov coefficients. It is, however, not necessary to perform the details 
of such an analysis because all the three spacetimes (Schwarzschild, de Sitter 
and Rindler) have virtually identical kinematical structure. In the case of 
Schwarzschild metric, it is well known that the thermal spectrum at late times 
arises because the modes which reach spatial infinity at late times propagate 
from near the event horizon at early times and undergo exponential redshift. 
The corresponding result occurs in all the three spacetimes (and a host of 
other spacetimes). 

Consider the propagation of a wave packet centered around a radial null ray in 
a spherically symmetric (or Rindler) spacetime which has the form in equation 
(162) or (175). The trajectory of the null ray which goes from the initial 
position at tin to a final position r at t is determined by the equation 




(186) 



where the • • • denotes terms arising from the transverse part containing dr^ 
(if any). Consider now a ray which was close to the horizon initially so that 

{f^in — I) I and propagates to a region far away from the horizon at late 
times. (In a black hole metric r ^ r^^ and the propagation will be outward 
directed; in the de Sitter metric we will have r -C Tin with rays propagating 
towards the origin. ) Since we have /(r) — > as r — > Z, the integral will be 
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dominated by a logarithmic singularity near the horizon and the regular term 
denoted by ■ ■ ■ will not contribute. [This can be verified directly from (162) or 
(175).] Then we get 



t-tin = ±[:r] / (1 + • • ■r'' dr^±[ — ]ln \f{rin)\ + const.(187) 





As the wave propagates away from the horizon its frequency will be red-shifted 
by the factor cu oc (l/y^^) so that 



9oo{rin) 
V ^00 (r) 



1/2 



J{r) 



(188) 



where K is an unimportant constant. It is obvious that the dominant be- 
haviour of u!{t) will be exponential for any null geodesic starting near the hori- 
zon and proceeding away since all the transverse factors will be sub-dominant 
to the diverging logarithmic singularity arising from the integral of (l//(r)) 
near the horizon. Since uj{t) oc exp[±(7t] and the phase 6{t) of the wave will 
be vary with time as 0{t) = J uj{t)dt oc exp[±g(t], the time dependence of the 
wave at late times will be 

'il){t) oc exp[i6'(t)] oc expi j w{t)dt oc expiQe^^* (189) 



where Q is some constant. An observer at a fixed r will see the wave to have 
the time dependence exp[i^^(t)] which, of course, is not monochromatic. If this 
wave is decomposed into different Fourier components with respect to t, then 
the amplitude at frequency v is given by the Fourier transform 

oo oo 

/(t/) = j dtij{t)e-"'' oc J e'^^'^-'^'dtoc J (itg-^^'^^-'^^^PM (190) 

— oo — oo 



Changing the variables from i to r by Qe^^^ — r, evaluating the integral by 
analytic continuation to Im r and taking the modulus one finds that the result 
is a thermal spectrum: 

1 Itt 

The standard expressions for the temperature are reproduced for Schwarzschild 
[g = (4M)^^), de Sitter {g = H) and Rindler spacetimes. This analysis stresses 
the fact that the origin of thermal spectrum lies in the Fourier transforming 
of an exponentially red-shifted spectrum. 
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But in de Sitter or Rindler spacetimes there is no natural notion of "energy 
source" analogous to the mass of the black-hole. The conventional view is to 
assume that: (1) In the case of black-holes, one considers the collapse sce- 
nario as "physical" and the natural quantum state is the Unruh vacuum. The 
notions of evaporation, entropy etc. then follow in a concrete manner. The 
eternal black-hole (and the Hartlc-Hawking vacuum state) is taken to be just 
a mathematical construct not realized in nature. (2) In the case of Rindler, 
one may like to think of a time-symmetric vacuum state as natural and treat 
the situation as one of thermal equilibrium. This forbids using quantum states 
with outgoing radiation which could make the Minkowski spacetime radiate 
energy - which seems unlikely. 

The real trouble arises for spacetimes which are asymptotically de Sitter. Does 
such a spacetime have temperature and entropy like a collapsing black-hole? 
Does it "evaporate" ? Everyone is comfortable with the idea of associating 
temperature with the de Sitter spacetime and most people seem to be willing 
to associate even an entropy. However, the idea of the cosmological constant 
changing due to evaporation of the de Sitter spacetime seems too radical. Un- 
fortunately, there is no clear mathematical reason for a dichotomous approach 
as regards a collapsing black-hole and an asymptotically de Sitter spacetime, 
since: (i) The temperature and entropy for these spacetimes arise in identical 
manner due to identical mathematical formalism. It will be surprising if one 
has entropy while the other does not. (ii) Just as collapsing black hole leads 
to an asymptotic event horizon, a universe which is dominated by cosmolog- 
ical constant at late times will also lead to a horizon. Just as we can mimic 
the time dependent effects in a collapsing black hole by a time asymmetric 
quantum state (say, Unruh vacuum), we can mimic the late time behaviour 
of an asymptotically de Sitter universe by a corresponding time asymmetric 
quantum state. Both these states will lead to stress tensor expectation values 
in which there will be a flux of radiation, (iii) The energy source for expan- 
sion at early times (say, matter or radiation) is irrelevant just as the collapse 
details are irrelevant in the case of a black-hole. 

If one treats the de Sitter horizon as a 'photosphere' with temperature T — 
{H/27r) and area Ah = AnH^"^, then the radiative luminosity will be {dE/ dt) oc 
T^Ah (X H^.li we take E = {l/2)H-\ this will lead to a decay law [334] for 
the cosmological constant of the form: 



k{L'pK)Uk,{t-U)) 



-2/3 



OC (Llt)-'/^ (192) 



where A; is a numerical constant and the second proportionality is for t — * oo. It 
is interesting that this naive model leads to a late time cosmological constant 
which is independent of the initial value (Aj). Unfortunately, its value is still 
far too large. These issues are not analyzed in adequate detail in the literature 
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and might have important imphcations for the cosmological constant problem. 



11 Cosmological constant and the string theory 

A relativistic point particle is a zero dimensional object; the world line of 
such a particle, describing its time evolution, will be one dimensional and the 
standard quantum field theory (like QED) uses real and virtual world lines of 
particles in its description. In contrast, a string (at a given moment of time) 
will be described by an one dimensional entity and its time evolution will be 
a two dimensional world surface called the world sheet. The basic formalism 
of string theory — considered to be a possible candidate for a model for 
quantum gravity — uses a two dimensional world sheet rather than the one 
dimensional world line of a particle to describe fundamental physics. Since the 
point particle has been replaced by a more extended structure, string theory 
can be made into a finite theory and, in general, the excitations of the string 
can manifest as low energy particles. This provides a hope for describing both 
gauge theories and gravity in a unified manner. (For a text book description 
of string theory, see [41,42]; for a more popular description, see [335,336,337]). 

It was realized fairly early on that string theory can be consistently formulated 
only in 10 dimension and it is necessary to arrange matters so that six of 
these dimensions are compact (and very small) while the other four — which 
represents the spacetime — are presumably large and non compact. There is 
no fundamental understanding of how this comes about; but the details of the 
four dimensional theory depends on the way in which six extra dimensions are 
compactified. The simplest example corresponds to a situation in which the six 
dimensional geometry is what is known as calahi-yau manifold [338,339,340] 
and the four dimensions exhibit N — 1 supersymmetry. The current paradigm, 
however, considers different ten dimensional theories as weakly coupled limits 
of a single theory and not as inequivalent theories. Depending on the choice 
of parameters in the description, one can move from one theory to other. In 
particular, as the parameters are changed, one can make a transition from 
weakly coupled limit of one theory to the strongly coupled limit of another. 
These strong-weak coupling dualities play an important role in the current 
paradigm of string theories though explicit demonstration of dualities exists 
only for limited number of cases [341,342,343,344]. 

The role of cosmological constant in string theories came into the forefront 
when it was realized that there exists a peculiar equivalence between a class 
of theories containing gravity and pure gauge theories. One example of such 
a duality [345] arises as follows: A particular kind of string theory in ten 
dimension (called type II B string theory) can be compactified with five of the 
dimensions wrapped up as 5-sphere (5^) and the other five dimensions taken to 
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describe a 5- dimensional anti de Sitter spacetime with negative cosmological 
constant (AdS^). The whole manifold will then be x AdS^, with the metric 
on the AdS sector given by 

ds'^ ^dr^ + e'^'~{r]i,^dx^dx'') /x, i/ = 1, 2, 3, 4. (193) 

This string theory has an exact equivalence with the 4-dimensional N — A 
supersymmetric Yang-Mills theory. It was known for a long time that the 
latter theory is conformally invariant; the large symmetry group of the AdS^ 
matches precisely with the invariance group of Yang-Mills theory. The limit 
r — > cxD is considered to be the boundary of AdS space on which the dual 
field theory is defined. This allows one to obtain a map from the string theory 
states to the field which lives on the boundary. It must be stressed that it is 
hard to prove directly the equivalence between type II B AdS^ x string 
theory and the four dimensional Yang-Mills theory especially since we do not 
have a non perturbative description of the former. In this sense the Yang- 
Mills theory actually provides a definition of the non perturbative type II B 
AdS^ X string theory. It is, however, possible to verify the correspondence 
by restricting to low energies on the string theory side. 

If gravity behaves as a local field theory, then the entropy in a compact region 
of volume will scale as S oc while indications from the physics of the 
horizons is that it should scale as oc R^. One can provide a consistent 
picture if gravity in D— dimensions is equivalent to a field theory in D — 1 
dimension with the entropy of the field theory scaling as the volume of the 
D — 1 dimensional space which, of course, is the same as the area in the original 
dimensional space. This is achieved in a limited sense in the above model. 

The AdS spacetime has a negative cosmological constant while the standard 
de Sitter spacetime has a positive cosmological constant. This change of sign 
is crucial and the asymptotic structure of these theories are quite different. 
We do not, however, know of any solution to string theory which contains 
de Sitter spacetime or even any solution to standard Einstein's equation with 
a positive cosmological constant. There are, in fact, some no-go theorems 
which state that such solutions cannot exist [346,347,348]. This, however, is 
not a serious concern since the no-go theorems assume certain positive energy 
conditions which are indeed violated in string theory. 

If de Sitter solutions of the string theory exists, then it would be interesting 
to ask whether they would admit a dual field theory description as in the case 
of anti de Sitter space. Some preliminary results indicate that if such a duality 
exists, then it would be with respect to a rather peculiar type of conformal 
field theories [349,350,351]. The situation at present is reasonably open. 

There is another indirect implication of the string theory paradigm for the 
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cosmological constant problem. The detailed vacuum structure in string the- 
ory is at present quite unknown and the prehminary indications are that it 
can be fairly complicated. Many believe that the ultimate theory may not 
lead to a unique vacuum state but instead could lead to a set of degenerate 
vacua. The properties of physical theories built out of these vacua could be 
different and it may be necessary to invoke some additional criterion to se- 
lect one vacuum out of many as the ground state of the observed universe. 
Very little is known about this issue [352] but advocates of anthropic principle 
sometimes use the possibility multiple degenerate vacua as a justification for 
anthropic paradigm. While this is not the only possibility, it must be stressed 
that the existence of degenerate vacua introduces an additional feature as re- 
gards the cosmological constant [353]. The problem arises from the fact that 
quantum theory allows tunneling between the degenerate vacua and makes 
the actual ground state a superposition of the degenerate vacua. There will 
be an energy difference between: (i) the degenerate vacua and (ii) the vacuum 
state obtained by including the effects of tunneling. While the fundamental 
theory may provide some handle on the cosmological constant corresponding 
to the degenerate vacua, the observed vacuum energy could correspond to the 
real vacuum which incorporates the effect of tunneling. In that case it is the 
dynamics of tunneling which will determine the ground state energy and the 
cosmological constant. 
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